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ON THE INFLUENCE OF THE DEMAGNETIZING FIELD ON 
DOMAIN STRUCTURE 


Bv J. KociŃski 
(Received February 4, 1958) 


It has been found that the demagnetizing field acting near the end of a crystal rod 
with mean magnetization parallel to the long axis of the rod, exerts a great influence on the 
domain structure of the crystal. The result of the approximate calculation shows, that this 
effect explains partly the discrepancy found between esxperiment and Néel's theory of 
the domain structure of such a crystal. It may be that the discrepancy that remains, de- 
pends on the approximations made in the calculation. The formula obtained for the double 
domain width d goes over into Néel's formula for a needle-shaped crystal. 


1. Introduction 


Néel (1944) found the domain structure of an iron crystal (Fig. 1), and investi- 
gated the change of this structure in an external magnetic field. But the predicted 
values of the double domain width d are nearly two times smaller, except in very 
weak fields, than those measured experimentally (Bates et alii 1956a). Іп Néel's 
treatment the crystal has the form of an rectangular rod whose length is great in com- 
parison with the two other dimensions. But this condition has not been fulfilled in 


Fig. 1. Domain structure of an iron crystal proposed by Néel. The external magnetic field H is parallel 
to the direction [110]. 


; (283) 
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the experiments. Therefore the demagnetizing field acting near the ends of the rod — 
which is neglected in Néel's consideration — may exert an influence on the domain 
structure and cause the discrepancy which has been found between experiment and 


theory. 


2. Qualitative consideration of the influence of the demagnetizing field 


During the measurement of the double domain width d, the crystal is placed in 
an uniform external magnetic field parallel to the long axis of the crystal rod. Then d 
is found as a function of the effective field (external minus demagnetizing). In the 
middle part of the crystal the demagnetizing field is nearly cancelled ош by the 
external field, but near the ends the domains are under influence of the demagnetizing 
field. Its direction is opposite to the direc- 
tion of the external field (Fig. 2). When 
the external field is increasing, the magne- 
tization vectors in the main domains, in the 
whole crystal, move towards the direc- 
tion of this field ; but this leads to an increa- 
se of the magnetic energy density connected 
with the demagnetizing field in the domains 
near the ends. Therefore the magnetization 
vectors will not move so much towards the 
direction of the external field as they would 
have done in the absence of the demag- 
netizing field. This will happen in the whole 
Fig. 2. Domain structure and magnetic field crystal, not only near the ends, since 
near the end of a crystal. otherwise magnetic poles should appear 
between the domains in the middle part and 
at the ends of the crystal. This hindering of the movement of the magnetization vec- 
: tors is equivalent to an increase of the width of the main domains, in comparison with 
the theoretical value, calculated with the neglect of the demagnetizing field. 


3. The energy of the domains at the ends of a crystal rod 


The inhomogenous demagnetizing field causes a space density of magnetic 
charge. We will neglect this space density and assume the domains to be magnetized 
uniformly. It has been shown by Williams, Bozorth and Shockley (1949) that a small 
space density of magnetic charge acts in the direction of decreasing the magnetic 
energy. This will be true also for large space density in the case considered. But 
since the final approximate result does not depend on the exact value of the magnetic 
"energy, the space density may be neglected. Then the demagnetizing field Н, is 
simply the field of a surface magnetic charge of uniform density I, cos O. It may be 
calculated as shown by Ferraro (1954 p. 150) (Fig. 3). 
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In the energy formulae the excess energy density of the closure domains Q, P 
(Fig. 1) over the main domains is important. With the components of Н. as in fig. З 
it may be seen that only the z component of Н. is important. We designate the depth 
at which H, is acting by д, (Fig. 2). The value of 6 will be found later. We will now 
consider the situation in the range of д. The intensity of the z component H4 of the 
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Fig. 3. Intensity of demagnetizing field along the lines A, B, C, D, for surface distribution of magnetic 
charge оп 51, Sg, with uniform density I, cos © = 1580 - 0.77 gauss. 


demagnetizing field decreases in the direction away the centre of a cross section of 
the crystal towards the edges (Fig. 3). Therefore we shall divide the crystal into three 
parts (Fig. 4). The effective field in the central part has the direction of the demagnetiz- 
ing field, in the two outer parts it has the direction of the external field. These later 
parts increase with the external field. 

Now we shall write down the expressions 
for energy analogous to those of Néel. In the 
central part (Fig. 4), we have the following 
situation: In the main domains the magnetic 
energy density is equal to: 


Е, = 3 L Ha cos O; (1) 


Нл stands for the intensity of the average 
demagnetizing field acting in the central part p ig. 4. Directions of effective magnetic field 


along the distance ô.. near the end of a crystal 
# 
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The magnetocrystalline energy density is: 
DNE. 
Ру == T К, (2 cos? 0 —1)*; (2) 
K, is the anisotropy constant, © is the angle between the direction of the external 


field and the magnetization vector, cf. Fig. 2. 
In the Q domains: 


Е„ = 0, Е, = 0. (3) 
In the P domains: 
1 1 1—£ © 

Fm = — < Is На созл = > 1,4—- dotz, (4) 

I Ө. m : ү я 
where Hj — 4 zj d cot PES the average demagnetizing field acting in the 

1 dlc 0 : 

domain. It depends on the height LS d cot 57 of the domain (Fig. 5), since the 


Fig. 5 Shape of the closure domains Q, P. 


intensity of the field given by the surface charge increases when we move from the 
edge to the centre of a cross section of the crystal. The parameter 4 determines the 
intensity of the demagnetizing field. 

The magnetocrystalline energy density is: 


=> К. (5) 
Thus the excess energy density in the Q domains is equal to: 


AFg = — E 1, На cos © + E К; (2 cos? O — y (6) 


we put: z I, На cos © = 0, 3 К, (2 cos? 9—1)? = Е, (7) 
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and in the P domains: 


1 1—5 OT 
ДЕр = = 1,4 оо, (8) 
" 1 1-4 9 1 
where we put: 3 І, А зом d cot ie P; ie SZR. (9) 


If the height of the central part is taken to be h(1—e) (Fig. 4), the sum of the 
excess energies in (© and P domains is: 


Е, = h(1 — s) B 62 d? tan © (CQ—F) + x (1— &)? d? cot E (Pat K-0- | ; 
| (10) 


1 
where A é?d? tan O, 5 (1 — £)? d? cot > are the volumes per unit height of 


the Q and P domains, cf. Fig. 5. 
In the outer parts we have the following situation. In the main domains the 
energy density in the effective field is: 


Ер = —I, H cos Ө (11) 


H, determines the average intensity of the effective field (external minus demagnetiz- 
ing) acting in the domains. 
The magnetocrystalline energy density is: 


Фу = i К; (2 соз? © — 1)?, (12) 

In the О domains: 
Ең == 0, F, = 0. (13) 

In the P domains: 


(again H, is the average intensity of the effective field), 
K,. (15) 


We put: І, Н, созӨ = а, ILH, = В. If h 2 designates the height of an outer 


part (Fig. 4), the sum of the excess energies in Q and P domains is equal to: 


F,— he E Pd tan 0 (a—F)+ 4 Й — E)? d? cot = бы ая n.| . (16) 


z 
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нн он‏ کے 
The whole excess energy in Q and P domains is equal to:‏ 


F = F, + Б = Фа O [fea В) — (1— 2) Q] + 


-la-see $ [0 — ®(Р — 0) + (ea— F)— (е8 – ЮП (17) 


per unit height. 

Since P contains £, this expression depends on the third power of &. We need the 
value of & for which it obtains minimum value. Since Ғ contains four unknown рага- 
meters only an approximate value of £,;, will be found. The condition: 


dF 
че 


leads to the equation: 
B (1— =) „| ё — n (l—e)y + c + cot а (ва — F) — 
0-90-90 (её — к) | +] жа—9я+ 9 te -0-0-30- 


mic be ку] -- 0 (8) | 


vith = I, Ad cot? (2) (19) 


When we put the discriminant of (18) equal to 0 the calculated value of Є, corre- 
sponds to the turning point of the curve Е(2), equation (17). It amounts to: 


cot — 
2 


ё, НОВ. (20) 
tan Ө + cot — 
2 
The values of $min and £,,, corresponding to minimum and maximum of F(£) are to 
be found at equal distances to the left and right of the turning point. They satisfy 
the condition 0< 6 < 1. This may be shown if one considers the expression (17) 
for F, where: Q > 0, F > 0, P > 0, and га < Q, and if one takes into account the 
fact that: 0,707 < & < 1 from (20). 
To get approximate results we shall assume that Є пах lies half way between &, 
and 1. We thus have for ć,,, the value: р 
1 2 cot Б — tan O 
(ep JAK WS LL Ro C YU (21) 
2 (cot э + tan ө) 
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The difference between this 6,,, and Є, amounts to 17% in a mean effective field of 
90 oersteds in the crystal and it tends rapidly to zero with an increasing field. There- 
fore in order to calculate the parameters Q and 7, we will make use of the condition 
that the discriminant of (18) is 0, which leads to the following equation: 


1 91 і 
Q2 |ia-e (une | со $) А 0)1 — =) [a (no + cor} AE 


+ (1— e) tan O — («В — K) cot З (во + cot 2) | de 


с с 2 
+ з[ва- 0 (шө +$) cot (В | -- 


+4 (l 95 (ca — F) tan © = 0 (22) 


Q contains the parameter H,,, which determines the intensity of the average demagnet- 
izing field. This is uniquely determined and therefore only one root of the equation 
(22) for Q has physical sense. To simplify the calculation we will put the discriminant 
of (22) equal to zero and will find by means of this condition the approximate value 
of 7 containing the parameter A. We obtain for О and 7): 


(ea — Е) (“зе + cot 3 + (ЕВ — К) "ве. 
B 5 2 
AE AW Е a ; (23) 


(1 — є) (te © + cot] 


8 С © + cot > cot E (eB — К) 


ар ТОЕЛ , > 
By means of 7 we obtain for P the expression: 
| 2 (sB — К) > 
EET: (52) 


Тһе value (20) for £, has been calculated by means of these expressions. Since, accord- 
ing to (19), 7 is positive, £8 — К must also be positive. We shall discuss this condition 


Jater. 


4. The Calculation of the double domain width d. 


We will find d by using the condition of minimum of the excess energy in the 
closure domains and the wall energy between the main domains. For the middle 
part of the crystal the expressions for the excess energy are these of Néel (1944). For the 
wall energy у we will take the expression given by Néel. When the range of the demag- 
netizing field is д we obtain the following energy expression: 

# 
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26 
32 sere | zii 


MEE + 5-8) (Р- О)--(ға- Е)-(еф--К)) 


per unit height, 


1 © 

Where: Sp = it xài: d? tan Ө, Sp = — (1 — %min)* d? cot — , 
4 4 2 
RZE MUT 9 

50 — 4 d : d? tan ©, Sp = = (1 mE Emin)” d? cot 9 » 


аге the volumes per unit height of the closure domains in the middle part and at the 
ends of the crystal, correspondingly, Fp and Fo are the excess energy densities in the 
closure domains of the middle part а to Néel, L, is the length and L the 
width of the crystal. The minimum of W is for: 


(26) 
h Wo = + О.Е W 1 F 
where о = т tan Q^ p^ po 5 
9 
tan O cot — : 
2 3 (1 — tan O 
4 (tun 6 + еи) 2 (ton Ө + cot 7) 


We will take the depth 6 from the field diagram (Fig. 6, 7). The zero point 
of the effective field in the middle part of the crystal rod is taken arbitrarily as the 
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Fig. 6, 7. Range д of influence of the demagnetizing field. Mean intensity of the effective field is measured 
from the lowest point of the upper curve — points Q, — and is equa lto about 50 oersteds on diagram 6, 
and 200 oersteds on diagram 7. 


lowest point ofthe upper curve. The value of 6 depends on the effective field 
approximately linearly according to the equation: д = c, Н, +c, with c, = 
= --0.065 cm/oe c, = 0.38cm found from the diagrams. 


6 
With growing crystal length L,, we obtain for d the formula of Néel, since 7- D dimi- 


nishes to 0. As we have said above, (ef da must be positive. Since in zero effective field 
Ер = 0, also К must be 0, But К = } Kj, as the anisotropy energy of the P domains. 


Table T 


3.7 


3.6 
0.80 3.3 0. 32 
0.85 2.9 0.26 
0.90 2.5 0.25 
0.95 2.0 0.24 


Tab. 1. The energy terms Wo, Wp and the double domain widths d, according to Néel's theory, and dą 
with the ends of the crystal taken into account. Calculated with K, = 2.8 x 105 erg . em", І, = 1580 
gauss for 3.896 Feb crystal length L, — 1.25 cm, and width L — 0.072 cm. 
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Therefore we must assume that instead of the P domains a different type of closure 
domains exists and the condition ef — K > 0 is fulfilled. The experimental evidence 
really reveals a closure domain structure somewhat different from that of P domains, 
with which another value of К is connected. Since the factor by which #8 — К is mul- 
tiplied is still smaller than 0.1 and г) is about 104 erg сш up to 105 erg спі З the 
second term in the denominator of (26) is substantially smaller than the first. When 
taken into account it would increase the value of d. But since (eB — К) may not be 
calculated precisely we will neglect this term while computing the values of d given in 


Tab. 1 and Fig. 8. 


120 4 4 
Experimental points from Bates and Hart (1956) 


—— — Neel’s theoretical curve 


== Theoretical curve obtained taking 
into account the ends of the crystal 


d (microns ) 


700 ' 200 300 
Effective field H (Qe) 


Fig.8. The change of the double domain width d in the effective field according to Tab. 1. 


5. Discussion of the result 


The discrepancy between experiment and theory that remains depends on several 
factors, apart from the approximations made in the calculation: 

1) The omitting of the second term in the denominator of (26); 2) the fact that 
instead of P domains other domains exist with presumably smaller excess energy and 
this causes a decrease of the first term in the denominator of (26) and further increase 
of d. This point is connected with a suggestion of Bates (1957); 3) the influence of the 
demagnetizing field on the wall energy. 

The discrepancy between Néel's theory and experiment should diminish with 
an increasing length of crystal with no change of its height and width. Since the 
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energy near the end is different from that in the middle part of the crystal, the values 
Xmin апа Emin for the same © are different. In small effective fields when the two types 
of closure domains are well visible, this difference should be detectable also for 
a short crystal, cf. Tab. 2. 


'Table II 


Tab. 2. The values of x,,;,, and Śmin determining the extension of closure domains in the middle part 
and near the end of the crystal. 


Bates and Hart (1956a) have found that the difference between experiment and 
Néel's theory decreases with the width L of the crystal. This should be expected 
because the intensity of the demagnetizing field decreases when Г, decreases. In 
another paper Bates and Hart (1956b) have found that the domain structure proposed 
by Néel is energetically less favourable than that given in their paper. Experiment 
speaks in favour of this new structure for small, but not for large crystal width Г. 
While taking into account the demagnetizing field it may be seen that in the case 
"of Néel's structure a smaller excess energy corresponds to the P domains (because 
the angle between the J,’s diminishes to zero), than to the Q domains (the angle be- 


ser л 
tween the J,’s is still equal to mak In the case of the structure of Bates and Hart, there 


are only closure domains analogous to Q, with rather large excess energy. This could 
be the cause of the fact that for a thick specimen Néel's structure is still energeti- 
cally preferable to the new one. 

I should like to thank Professor Szezepan Szezeniowski for helpful discussions 


and constant encouragement during the course of this work. 


КРАТКОЕ СОДЕРЖАНИЕ 


E. Коцинский, Влияние размагничивающего поля на областную структуру. 


Указано, что размягничивающее поле, которое действует вблизи концов кри- 
сталла намагниченного в среднем параллельно в его оси, в значительной степени 
влияет на его областную структуру. Приближенный расчет указывает, что это 
явление частично выясняет несогласие между опытными фактами и теорией 
Неэля изложенной ним для кристаллов железа. Быть может, то что согласие 
между теорией и фактами еще не является полным, связано с приближенным 
характером произведенных исчислений Получается формула для удваенной 
толщины d областей, которая переходит в формулу Неэля для иглообразных 
кристаллов. 
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MEASUREMENT OF ABSORPTION IN GRAPHITE OF ELECTRONS 
FROM MESON „ DECAY 


Ву W. Козсн 


Institute of Experimental Physics, Łódź University 


(Received November 25, 1957) 


Measurements were made of the absorption in graphite of electrons from the decay 
of д mesons. The results were compared with the absorption curves calculated by assum- 
ing a linear relation between the actual electron path and the energy, and relations that 
result from inclusion of the radiation losses. In the calculations the electron spectrum 
given by Michel (o = 0.55, W = 52.7 MeV) was assumed. Account was also taken of the 
increase in the path of the electrons as a result of scattering. Comparison of experimental 
and calculated results show that the mean total energy losses for electrons from the decay 
of д mesons is 1.92 MeV/g cm^?, whereas the mean energy losses due to ionization is 


(1.70 == 0.15) MeV/g cm-?. 


I. Introduction 


Ап attempt has been made in the present paper to determine the average total 
ionization losses for electrons from the decay of mesons u absorbed in graphite. 
The theory of the ionization loss by charged particles passing through matter 
was worked out by Bohr (1915), subsequently developed in the thirties by Móller 


(1932), Bethe (1932), and Bloch (1933). The result of theoretical calculations of the 
mean ionization losses 


-ELN | во(рав 5 (1) 


(where М is the number of ® electrons in 1 cc of absorber, с (Е) dE the cross section 
for collision with a transfer of energy between E and E + dE) is dependent on the 
` correct determination of the minimum, and maximum energy which can be transferred, 
and on the assumption of the correct cross section. Swann (1938) pointed out that 
in condensed media the effect of the mutual screening of the individual atoms should 
limit the relativistic increase of the energy losses. А simplified theory of this effect 
was given by Fermi (1940) and was later developed by Halpern and Hall (1949), 
Budini (1953) and Sternheimer (1954). The fluctuations in the energy losses during 


; (295) 
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travel through thin layers of absorber were considered by Williams (1931) and Landau 
(1944). Experimental work on the measurement of ionization losses has been carried 
out by various methods. Аз regards the measurement of absorption for the purpose 
not only of determining the energy losses but also the range of charged particles (this 
is of practical importance) the results of work done prior to 1952 have been described 
in an article by L. Katz (1952) and A. S. Penfold (1952). It was shown that the range 
of electrons given in Al on the basis of measurements by various authors for the 
energy interval of 1—16 MeV can be given by the formula: 


R (g/cm?) = 0.530 E, (MeV) — 0.106 (2) 


It was also found that the value of dE/dR is about 25% greater than the ionization losses 
predicted by theory, however, the reason for this difference was not established. 
The results of the work of Hereford and Swann (1950) and Beatti (1951) do not con- 
firm the difference between experiment and theory. An attempt could be made to 
resolve this disagreement on the basis of data obtained by other methods, but each 
of these methods has certain shortcomings that limit the accuracy. (А review of 
methods of measurement and an analysis of the results obtained up to 1956 have 
been given by Price (1956)). Measurements made by the proportional counter and ioniza- 
tion chamber methods show a relativistic increase in ionization, but it is only in the 
experiments of Parry (1953) and Eyeions (1955), carried out with pressure chambers, 
that the slope of the experimental curve for u mesons was found to be qualitatively in 
accordance with the theory of Fermi. In measurements made with Wilson chambers 
the cluster effect hinders the accurate determination of the mean ionization. The 
paper by Ghosh (1954) is one of the few in which the moments of the mesons have 
been determined with sufficient accuracy and an analysis made of the fluctuation 
effects. It shows that the proper ionization on the plateau is in accordance with the 
predictions of Budini's theory. The method of low-pressure counters contains errors 
due to knock-on electrons knocked out from the walls of the counters. There are lar- 
ge differences between the results of various authors. The results of Hereford (1948) 
. speak in favour of the Halper and Hall theory whereas the results of McClure (1953) 
lie on the Bethe curve. As for the photographic emulsion method, it can be seen on 
the basis of the comparison in Price's article, which includes experimental points of 
many authors, that there is а 15% difference between the ionization minimum and 
. the plateau. In view of the spread of the various results, however, and the fact that 
there exists a threshold energy below which a grain is not developed it is not possible 
to obtain an answer as to quantitative agreement or disagreement with the various va- 
riants of the theory. In measurements of energy losses in crystals and liquid scintilla- 
tors the observable quantity is the light output, which, in the first approximation, 
depends linearly on the ionization losses up to energies 4.I,;. (Chou 1952). Recently 
Bowen (1954) found а 10% increase in ionization for relativistic particles іп Ма]. 
The method of measuring low-energy losses in thin foils by means of a magnetic | 
analyzer was used by Goldwasser (1955) for investigation of the density effect (15.7 MeV 


Measurement of Absorption in Graphite of Electrons н 297 


electrons) and by Hudson (1957) who made measurements by means of electrons of an 
energy of 150 MeV. From the results it can be said that the probable losses and the 
spectrum of the electrons behind the absorption foil generally are in agreement with 
the predictions of the Landau theory with Sternheimer's correction. In view of the 
. differences between the measurements obtained and the sparse data concerning the 


energy —— > 40 the study of this question seemed of interest. The mean energy of 


oc 
the electrons emitted during и meson decays is nearly twice the threshold energy of 
the Fermi plateau; it would therefore seem that despite the general shortcomings of 
the measurements made by means of particles with a continuous energy spectrum, 
sufficiently accurate data could be obtained on the magnitude of the mean ionization 
losses, since, for most electrons from u meson decays, these losses should be the same 
and schould be constant. 


II. Apparatus 


The experimental set-up (Fig. l) consisted of a three-tray vertical telescope 
emitting a stream of u mesons. The С. M. counter trays registering the decay electrons 
are located symmetrically around the telescope. Over the telescope is a 1 cm iron 
plate and 7.5 cm of lead to eliminate the soft component of cosmic radiation. A block 
of graphite measuring 5 х 20: x 40 cm was mounted underneath the telescope. In 
passing throught the block, the low-energy mesons, after slowing down, stop and 


ГА 
| 
M 


Fig. l. The geometry of the vettical telescope trays, and trays registering the decay electrons. 
# 


\ 
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—— 


decay. The trays consisted of G. M. counters of the Maze type with a diameter of 
2.8 cm. The active length of the counters at trays I, A, B, C, D was 50 cm. They 
were energized from a rectifier with electronic stabilization. The telescope counter 
pulses after amplification and shaping were fed through the collectors to the coinci- 
dence circuit R, with a resolving time of 0.5 ү sec (Fig. 2). Triple coincidences of the 
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Fig. 2. Block diagram of the electronic circuit. WF — amplifying and pulse forming circuits, K Ex pulse 


collectors. К — coincidence circuits. М — numerators. L — delay line. D — separating diodes. 


telescope were registered by numerator N,. A detailed analysis of the operation of the — 


various individual components of circuits of this type has been made in the paper of | 


A. Zawadzki (1954). The narrow pulse arising at the output of R, was fed from the 
cathode follower through a delay line (1 sec) to a second coincidence circuit К». 


It entered the ,,broadening,, channel in which it was formed into a rectangular pulse : 


with a width of 5 u sec. If, during this pulse, R, received two pulses from trays 4 


— 
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or B and C or D, then there occurred a coincidence (I + II + II) + AB + CD, 
registered by the numerator №, which was interpreted as the stopping of the 7 
meson in the graphite block, and a decay with the emission of an electron which 
reaches the side trays. The delay in the pulse from the telescope was necessary in 
view of the great probability of registrations of this type being induced by the meson 
and the knock-on electron knocked out of the lead or the iron or possibly by a local 
shower. 


III. Control Measurements 


The maintenance of the constant efficiency of the circuit and the prevention 
of an equipment background was the fundamental tasks of this work in view of the 
small counting rate, of the order of 1 per hr without absorber. 

The following methods of apparatus control were applied: 

а) control of the shape of the pulses at various points of the electron circuit by 
the use of a pulse generator and an oscillograph; 

b) control of the number of pulses in the auxiliary statistics: measurement of 
the number of mesons passing through the telescope by means of coincidence cir- 
cuits R; and R, (1390/hr), measurement of the number of coincidences in statistics | 
(I + II + III) + AB + CD and Ш + AB + CD without line delay; the numbers 
of registrations were 40/h and 900/h with absorbers and in both cases they diminished 
as the absorber thickness was increased (table I); 

c) control of the counter condition. 

The control measurements, made twice per 24 hours, required about 4 of the 24 
hours of the operation of the circuit. The apparatus worked without interruption for 
200 to 300 hours. 

In order to obtain information. about the quality of the apparatus, mainly about 
its suitability for operation during the extended period of time for which it was ex- 
pected to function, measurements were made of the mean life of the mesons. This was 
reduced to the registration of the coincidences (I + II + Ш) + AB + CD (without 
absorbers) with several delays on the delay line. The background 0.027/h was deter- 
mined by measurements described further on in this paper; provisional measurements 
and a calculated estimation that took into account the resolving times of the coinci- 
dence circuits, however, showed prior to the beginning of these statistics, that it is 


greater than 0.1/h. Fig. З gives the experimental results and the relation N = № e * 
for т = 2.15 u sec. 

The results for delays less than 1.2 u sec show increasing deviation from the ex- 
pected exponential curve resulting from the curve, this is due to knock-on electrons and 
local showers. On the basis of the data obtained from these measurements it was found 
that the apparatus was suitable for registration of coincidences occurring with a fre- 
quency even of the order of 0.1/h. After completion of this measurement the resol- 


ving time was reduced for the coincidencesbe tween channels I and II of circuit. Ка to 
р | 
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0.3 u sec; this made it possible to eliminate the background with a minimum delay 
of 0.8 usec and hence allowed the raising of the recording rate for electrons. The 
duration of the rectangular pulse in channel III of К, was increased from 3.5 to 5.0 usec. 
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Fig. 3. Results of the measurement of the mean life of м mesons in graphite. The function 


N = № exp — normalized for t = 1.2: 10-5 sec to N = 0.700, is shown as a solid line. 


IV. Measurement of Background 


Six measurements of the background were taken during the intervals between the 
basic measurements. They consisted in the registration of the coincidences (I + П-- 
+ Ш) + AB + CD with a delay of 0.8 or 1.0 u sec in the pulse from the telescope 
without the graphite block for slowing down the mesons. In five of the cases there 
were no absorhers between the side trays and in one case the measurement was made 
with a 10 cm graphite absorber. А mean value of 0.079 + 0.016/h was determined 
with a confidence level of 0.22. With absorbers, however, only one coincidence was 
recorded in 85.75 hours. A statistical analysis by means of Pearson’s funcction y? 
showed that this result does not fit within the limits of the fluctuation of the pre- 
vious five measurements. The number of random coincidences for (I + II + IIT) + 
+ AB + CD was determined on the one hand by interpolating of the results of auxiliary 
measurements: I + AB + CD (9.8 + 1.2/h) and III + АВ + CD (0.67 + 0.18/h) with 
a delay of at least 1.0 sec in the pulses from the telescope, the pulse rate from trays I 
and III being known; and on the other hand it was calculated from the double- 
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-coincidence rate AB -- CD (7000/h), the resolving time of the circuit R, (5 u sec) 
and the pulses rate I + II + Ш (1390/h). Both estimates place the number of random 
coincidences at 0.027/h. On this basis И was possible to explain the difference in re- 
sults between the five measurements with absorbers and the sixth without absorbers. 
Even with the graphite block removed a small number of mesons decayed in the layers 
of wood supporting the counters and in the walls of the counters of trays C and D. 
This factor from the decay electrons of the order of 0.05 disappeared (last measure- 
ment) when 10 ста of graphite was placed between the trays. The estimations taking 
into account the geometry of the set-up, the thickness of the layers of wood and the 
thickness of the counter walls yielded a result of the same order as the above figure. 
Tt was further assumed that the number of random coincidences was 0.027/h and that 
the electronic circuits operated without apparatus background. 


V. Measurement of Absorption in Graphite 


Measurements мөге made in three runs which were denoted by Р, О, and R. 
The most important changes made in the apparatus, e. g. change of counters, change 
of anode voltage or valve polarization voltage, change of a large number of valves, 
change in the delays on the line, made each run different. In order to reduce the 
systematic errors the measurements were made “simultaneously“ for two thicknesses 
of absorbers, every 24 hours the absorber thickness was changed. The results of the 
basic and auxiliary measurements were averaged out within the run for the given 
absorber thicknesses (d = 1.67 g/cm?) as shown in table I. As mentioned above, 
auxiliary statistics (make it possible to evaluate the operation of the entire circuit 
during the measurements. The table presents mean values from the individual measure- 
ments made during the periodic controls fór auxiliary statistics (I + II + Ш) + 
-- AB + CD and III + AB + CD without delay. Registrations of the first type were 
produced mainly by knock-on electrons (triggering of trays А and B) and hence 
their number depends more on the absorber thickness than in the case of the registra- 
tion of the second type, which could be produced by single mesons. The results of 
run R are displaced downward parallelly to the results of the run P and Q for measure- 
ments in which the particle producing the coincidence had to pass through the three 
trays of the telescope, whereas in the case of the coincidences Ш + AB + CD the 
results of the three runs agree within the limits of error. The relation between the 
registrations and the voltage of the counter power supply was determined on the basis 
of an analysis of the remaining auxiliary measurements and also several times at va- 
rious periods of the operation of the experimental arrangement. It was found that 
for coincidences of varios types this displacement was the result of a rise in the thres- 
hold of the counter plateau in connection with the aging of the counters. The thres- 
hold voltage for registration with full output of coincidences III -- 4B -- CD was 
about 40 V lover than in the case of (I + II + III) - AB + CD and this change did 


not affect its results. The contribution from a possible change in the output of the 
Ў 
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Table I 
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Tab. 1. The results of measurements in runs P, С, В, R,, В, — the mean number of registrationsh of 

ji mesons measured by means of the coincidence circuit Ry or R,; III + АВ-- CD and (I + 

„+ II + Ш) + АВ + CD — the mean number of coincidences of the third tray or three trays of 

the telescope with the side trays without pulse delay. Counter trays — mean background of trays C, D, 

A and B counts per minute with a 1:64 scalar. Result — the numerator indicates the number of 
registrations, the denominator the counitng time in hours. 


electronic circuit could not exceed 3%. The measurements were relative and the 
above-mentioned displacement of the curves was parallel hence, treating the various 
individual runs as being independent their mutual normalization was performed after 
eliminating the background by the method of least deviation. Then, on the basis of 
the results of the various runs, the weighted means were determined for various 
thicknesses of graphite and corrections were introduced for absorption in the counter 
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walls and the thin layers of wood supporting the counters (equivalent of 0.72 cm of 
graphite). Finally, normalization was carried out by assuming the counting rate for 
electrons without absorber as 1000/h. The error indicated in Fig. 4 is the sum of the 
statistical errors without background and the expected background. 


Fig. 4. The results of measurements of the coincidences of the side trays as a function of the graphite 
absorber thickness. Crosses indicate the points calculated by setting а = 0.46 and а = 0.50. 


VI. Determination of the Shape of the Expected Absorption Curve on the Basis 
of Data on the Spectrum, Scatter and Electron Range 


In this work the range of the electrons cannot be determined directly on the 
basis of an experimental absorption curve owing to the untypical geometry of the 
recording circuit and the large dimensions of the graphite block which was the source 
of the electrons. In order to obtain a result which could be compared with experi- 
mental data, the expected number of registrations were calculated for various thick- 
nesses of absorber by assuming a linear relation between the real path of the electrons 
in the graphite and the initial energy: Z = aK). The interpretation of the measure- 
ments was therefore reduced to the determination of the value of a by comparison of 
the experimental absorption curve with the analogous curves calculated for several 
values of this parameter. The calculations were based on the following reasoning: 
each element of the graphite block is a source of electrons with a spectrum given by 
the formula of Michel (1950). 

AE? 4 
var = 48 |з(7—Е) +20 (+ 5-я | dE (3) 
(where W is the maximum energy of the electrons, 0 a parameter of the spectrum 
shape). If an absorber is расей between the side trays then the inner side trays will | 


304 W. Kusch 


be reached only by electrons with an initial energy of E >> Ek, whose range is greater 
than the actual path in the graphite. The number of these electrons is: 


W 
[Е (КУПЕ = eee со Z EE (4) 


Ek 
The actual path of the electron depends on the thickness of the graphite layer. D. 
By setting Z = аЁ we have 


w 
[ ® dE = L= (2) &.- tpa (5) 
Ер 


In order to carry out the integration over all the elements of the graphite block (with 
coordinates x, y, z) it is necessary to take into account the fact that the number of 
mesons An* (xyz) dxdy dz decaying in the various elements is variable and the number 
of electrons counted is dependent on the solid angle Q (хуг) of a given element on the 
side trays 


М (D, a) = Јао һ — a Kı — (2) x] dx dy dz (6) 


The integration is carried out over the volume c of the block included in the stream 
of u mesons. 

I. In order to determine the values of An* (xyz) dx dy dz the following calcu- 
lations were made: 

а) An(xyz) dx dz the number of mesons м passing in a unit of time through 
any horizontal element of surface contained in the graphite block was determined. 

b) The mean path L of the и mesons in the graphite block to any element 
dx dy dz of the graphite block was calculated. In these calculations the directional 
distribution cos?g was taken into account. 

c) The distribution of the mesons stopped i in the individual elements was deter- 
mined. 

Assuming that for mesons of energy up to 120 MeV (the energy which the mesons 
lose in the graphite block) the relation ес the range and the energy сап be ap- 
'proximated by the function 


R=aE' MO 


where a is a constant, and that the low-energy part of the spectrum (0—120 MeV) об: 
the mesons at sea level can be approximated by the function 


N(E) dE = bE" dE (8) 
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where b is a constant, we obtain the distribution of the stopped mesons: 


N(R) dR = ск” dR (9) 
where c is a constant. 

Since the mean path of the mesons in the graphite is proportional (to an accuracy 
of the order of 1%) to the vertical path L of the mesons in the graphite, then the 
distribution of the mesons stopped along the vertical path can be represented by this 
same function. Hence the number of electrons emitted from the given element is 
expressed by 


An* (xyz) dx dy dz = KAn(xyz) y^ dx dy dz (10) 


where K is a constant 

II. The solid angle from the individual elements on the outer side tray was 
determined. Under the assumption of an isotropic distribution of the emitted elec- 
trons the number of electrons registered by the side trays should be proportional to 
An* Q. 

III. The actual path of an electron Z (xyz, D, E) was determined with the following 
assumptions. 

a) the mean path S(xyz, D) of an electron from a given element to the side trays 
was calculated by assuming an isotropic distribution. 

b) An estimate was made of the average increase in the path as a result of the 
scattering of the electrons. The basis for the calculations was provided by the work of 
Yang (1951) Actual Path of Electron in Foils in which the author gave.the value 
of the actual path Т = t + A as a function of the foil thickness t (in cascade units) 
and the energy of the electron 4 


_ № "an 
El uż where w = з (Ein Ме У). (11) 


A generalization was made for thick absorbers by applying Yang в formula to 
imaginary thin layers of graphite and taking into account the linear energy losses 
Е = E,—aT. These calculations were made once with the assumption that the 
electron loses energy in steps after passing through the foil and again with the assump- 
tion that it loses energy before entering the foil; this made it possible to find an 
approximation formula expressing the relation between the actual path and the thick- 
ness of the absorber for an electron that loses all of its initial energy E in the absorber 


ERSTE Ы” 
ано ра 19 
t Е ча 


Аз a result of the adoption of this approximation we obtain a relationship between 
the initial energy of the electron and the maximum absorber thickness through which 


zaj GI 


it can travel: 
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In further calculations the actual path of the electron was assumed to be equal to the 
mean path to the side trays determined from the element from which it was emitted, 
a correction being made for the increase in the path due to scattering. Hence Eq. (6) 
can be written in the form | 


якору DZE 3 
Noa- [ао [50 pis) а 
ie 4 
" A B E «yh ж See) Г dy dz (14) 


The evaluation of the above integral was done numerically. The graphite block from 
which the stopped mesons emitted electrons was divided into 1000 volume elements 


measuring 2 x 2 x 1 ст (where the dimension perpendicular to the side trays was 
1 em), the integrand was determined for the centres of the various elements, the sum 
being obtained on an adding machine. For а = 0.46 the points D = 0,6,8 and 10 cm 
were calculated. For the points determined in this manner a correction was made for 
the adoption of a mean path for all the electrons. This correction was determined by 
comparing the absorption curves calculated for one volume element under the assump- 
tion of an isotropic distribution, an angular distribution 1 + a cos $ (a — —0.20) 
and all the electrons travelling the same path, i. e. the mean paths. Experimental 
results and segments of the absorption curves are shown in Fig. 4. 

If calculations are made of the absorption curves for one volume element with 
а = 0.48, 0.50, 0.52, 0.54 then the ratios of the distances of the points of the indivi- 
dual curves will be the same as in the case of curves calculated for these same values 
of the. parameter for all volume elements. On this basis it was found that the best 
agreement between the experimental results and calculations is obtained for а = 0.52. 
This value is the reciprocal of the total losses of an electron over the actual path 
E В; 


ar RE “я 
з dz 1.92 MeV/g ст”. 


VII. Comparison of Experimental and Theoretical Results 


In order to compare the experimental results with the results of the theory of 
ionization losses it was necessary to determine the contribution of radiation losses. 
On the basis of the formula 


2E 1 -i 
іп س‎ 2 < E <> 137 2.7 
= _4DŻ. ТЖ "mad 15 |" зів 322 


CCL BOR Kya % 2 (15) 
à тів Да E > 137 m, c? Z 


Gol 


(where N is the number of atoms іп 1 ст, the atomic number of the stopping medium, 
го the classical radius of the electron) it is found that in the case of graphite, for 
energies above 40 MeV, the losses are linearly dependent on the energy and, for 
lower energies the bracketed factor depends very weakly on then energy. Therefore in 
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AM ae dE 
approximation we can take — | — = В. E(B = 1.74 х 10-2 cm?/g). Since for 
dO | ыш 
relativistic electrons the ionization losses can, with little error, be taken to be cons- 


dE dE 
tant — (зг). = A, then — 26]. = А + B: E, which leads to the relation 
between the path travelled and the energy: 


1 В 
z= pres = (16) 
This formula was applied in calculating the absorption curve; in the calculations, the 
increase in the path due to scattering was taken into account for three elements of 
volume (lying on an axis perpendicular to the side tray — coordinates z — 2, 0, 
—2 cm) and with a fixed parameter B for A = 1.54, 1.67, 1.82. As may be seen 


from Fig. 5 the experimental curve is contained between curves of a parameter 


1 


; | В г 
Fig. 5. Absorption curves calculated by assuming Z = T 1+ P в) for В = 1.74 : 10 * 


cm?/g, А = 1.54, 1.67, 1.82. Crosses mark the experimental curve. 


А = 1.67 and 1.82. On the basis of the above estimation the mean value of the ioniza- 
tion losses in graphite for electrons from the decay of u mesons is А = 1.70 F 
F 0.15 MeV/g ст”?. Since the mean radiation losses are 13% of the ionization losses, 
an estimate was made of the possibility that Bremsstrahlung quanta produce pairs 
and Compton electrons that might reach the side trays despite the use of an absorber 
thickness equal to the maximum range of the primary electrons. The probability of 
such a phenomenon for electrons of an initial energy of 30, 40 and 52 MeV is 2 х 10-3, 
3.6 x 10-3 and 6.9 x 10-3, respectively. This effect could have an influence on the 


results of the measurements. 
; 
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VIII. Conclusion 


Comparison of the results of the calculations made under the assumption that 
there exists a linear relation between the actual path of the electron and its energy 
shows that the best agreement between the calculated and experimental curves is 
obtained for a = 0.52, which corresponds to the mean total losses 1.92 MeV/g сп”. 
Computational errors resulting from the approximation, the assumption of an iso- 
tropic distribution of the emitted electrons and the uncertainty in the parameter o 
can change the value of а by not more than +2%. The assumption that the relation 
between the initial energy and the path travelled is a linear one justified on the one 
hand by the experimental results (Eq. 2) and on the other hand by considerations 
concerning the range with the radiation losses being taken into account; Eq. (16) can 
be approximated by a linear relation with a mean error of 2.5%. The determination 
of the radiation losses enabled the estimation of mean ionization losses. The result 
1.70 + 0.15 MeV/g cm”? is in agreement, within the limits of error, with the losses 
predicted on the basis of the Halpern-Hall formula. 

It seems that the difference between the experimental and theoretical results, as 
noted in the Katz and Penfold article comes mainly from the failure to take into 
account, for lower energies, the effect of the increase in the actual path of the elec- 
trons as a result of scattering and, in the case of relativistic energies, the failure to take 
radiation losses into account. 

I should like to express my thanks to Prof. M. Miesowicz, who directed this work, 
for valuable discussions and advice and to Docent A. Zawadzki for his advice and 
assistance during the work. I also wish to thank the Department Head, Prof. L. Na- 
tanson, for his interest and encouragement. 


КРАТКОЕ СОДЕРЖАНИЕ 


В. Куш, Измеренйе абсорбции в графите электронов с распада и-мезонов 


Произведено измерение абсорбции в графите электронов с распада ц-мезонов 
а результаты сравнено с кривыми абсорбции подсчитанными для. предположения 
линейной зависимости между действительным путем электрона и энергией, а так- 
же зависимости с учетом потерь энергии излучения. В расчётах принято спектр 
электронов по Мишелу (о--0,55; W=52,7MeV). Учтено также ‚Удлинение пути 
электронов по поводу рассеяния. Сравнение экспериментальных ‘и:подсчитанных 
результатов позволило заключить, что средние потери энергии для электронов 
с распада мезонов становят 1,92МеУ/гсм-?, а средние потери энергии через иони- 
`зацию (1,700,15)MeV/rcM-?2. 
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DISCUSSION OF THE TECHNIQES EMPLOYED IN THE MEASUREMENT OF THE ACTIVITY 
OF NATURAL CARBON 


(Received January 4, 1958) 


Comparisons are made of the techniques of measuring the radioactivity of natural 
carbon a) using solid carbon and screen wall counters b) using gas samples. The gas 
sample technique is advantageous from the point of view of sensitivity, use of material, 
change of samples, and less contamination. One difficulty with this technique when pro- 
portional counters are employed, is the complicated electronic apparatus and above all 
the lack of control of the background level during the measurement. Statistical methods, 
of controlling the background by employing the fluctuation of the intensity of the cosmic 
radiation in the counter and the shielding are discussed. 


1. Introduction 


The generation of the isotope !*C in the atmosfphere was predicted by Libby 
(1946). Work on the radioactivity of biomethane and methane that accompany petro- 
leum demonstrated the correctness of Libby’s hypothesis (Anderson et al. 1947, 
1947a, 1947b). Attention was drawn to the possibility of using the 14C generated in the 
atmosphere for chronometric purposes (Libby, Anderson, Arnold 1949; Arnold, 
Libby 1949). The technique measuring the concentration of !4C'in natural carbon 
by the use of screen wall counters (Libby & Lee 1939) and carbon in the solid state 
was developed by Libby and his co-workers. (Anderson & Libby 1951; Anderson, 
Arnold, Libby 1951a). By means of this technique a number of control measurements 
were made (Arnold, Libby 1949) and the age of samples of the late pleistocene period 
was determined. 

The suitability of the 14С method for the determination of the age of organic 
remains, the scope of the method, the possible sources of error, the dependability 
of the results, etc. һауе already been discussed (Flint & Deavey 1951, Flint 1951, 

Johnson et al. 1951, Godwin 1951, Anderson, Levi 1952, Mościcki 1954, and others). 
| ke 
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Independently of the work of Libby et al. the author, basing himself on the work 
of Brown and Miller (1947), proposed late in 1948 a technique employing samples - 
in the gaseous state. (Mościcki 1949—1950). An apparatus based on this method was 
assembled and put into operation in 1950—1951 at Poznań. Berendsen and de Vries 
(1952) developed a technique using pure CO; in a proportional counter. Barker (1953) 
developed a technique for the synthesis of C,H, from CO, obtained from the combus- 
tion of the samples, and Crathorn used C,H, in a proportional counter for measuring 
the 14C concentration (Crathorn 1953). Faltings (1952) considered the possibility of 
employing a technique in which methane is used. On the other hand, an essential 
improvement was made to the elemental carbon technique in the form of a double 
divided cathode which permitted the elimination of пш pulses during the 
measurement (Double S. W. C.). 

In recent years major progress has been made in the perfection of the technique 
for measuring the activity of samples and in the improvement of the methods employ- 
ing gas samples (е. g. Braunon, Taggart Jr. & Williams 1955; Burke Jr & Meinschein 
1955; Houtermans & Oeschger 1955). 

According to a report at the Symposium in Copenhagen, 1954, there were in 
Europe, in 1954, eleven absolute geochronology laboratories.! 

Four laboratories employed the technique using solid state carbon, the others 
the gas sample technique. Comments concerning the advantages of the solid-state 
carbon technique can be found in a paper by Anderson and Levi (Anderson & Levi 
1952). De Vries and Barendsen, Crathorn (l.c) and Mościcki (1953) among other 
authors, drew attention to the advantages of the gas counter technique over the solid- 
-state carbon technique. | 


2. General Remarks About the Advantages and Shortcomings of the Methods? 


The advantages and shortcomings of the method can be decided by: 

1. The dependability of operation of the 14С decay-electron detector. 

2. The sensitivity of the method and hence the minimum mass necessary for 
making a measurement with a certain statistical accuracy. 

3. The chemical processing, in particular: 

a) The ease with which the chemical compound required for the measurements 
may be obtained and the samples changed. 

b) Certainty of avoiding, during chemical processing, contamination with mate- 
rials of a radioactivity other than the test sample. 


w 


1 The Poznań apparatus was disassembled in connection with the author's transfer to the Nicholas 
Copernicus University in Toruń (1952—1954) and the Gdańsk Politechnical College (1954). 

? With the exception of the statistical distribution law of the ratio АМЫ AL; , the remarks in this 
part were presented in a paper read by the author at a Colloquium of the Geology Committee of the 
Polish Academy of Sciences in Torun (March 23, 1953). 
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с) Standardization of the counter gas in the sense of it being possible to reproduce 
the effective mass of the carbon introduced into the counter from measurement to 
measurement, the degree of purification from possible impurities, etc. 

4. Control of the counter efficiency and the stability of operation of the appa- 
ratus. Ease and accuracy of elimination of possible fluctuations in the background 
during the measurements. 

5. Simplicity of the electrical equipment. 

6. Size of the apparatus capable of making a given measurement. with a given 
statistical error. 

In view of the basic similarity of the gas counter techniques we shall centre our 
attention on the comparison of these techniques with the technique using carbon in 
the solid state. 

Ad. 1. The dependability of the detector operation is, to a great degree, a func- 
tion of the standardization of the conditions in which the apparatus operates; generally 
speaking, these conditions can be maintained the same for all techniques. Information 
published hitherto concerning the hehaviour of the detectors employed do not yet 
make it possible to draw final conclusions as to possible irreversible changes occurring 
in the counter which could affect the level of the background pulses (the number of 
delayed pulses, etc.) and as to irreversible changes in the counter efficiency 
from measurement to measurement owing to ageing or other possible effects that 
distort the results of the measurement. For this reason we must, at least for the pre- 
sent, recognize all techniques as being of equal value as far as dependability of detectors 
is concerned. It must be borne in mind, however, that the properties of the argon- 
-alcohol counters used in the solid carbon technique are well known and no unpleasant 
surprises are to be expected in this case. However the properties of proportional 
counters filled with СН,, CO,, C,H, and particularly the counter filled with CO, + 
--CS, working in the С. M. range are rather less known and for this reason the equating 
of the reliability of operation with that of argon-alcohol counters is, in a way, confi- 
dence on credit. This applies in particular to counters filled with CO, + CS. The 
papers of Miller and Miller et al. (1950) showed that these counters do not give the full 
efficiency ; their efficiency was estimated to be 85%. Direct measurements of the 

efficiency made by Hawkings, however, indicated nearly 100% efficiency, and the 
work of Mann and Parkinson (1949) did not show a relatively low efficiency for the 
counter. Work by the author in Poznan and Torun in the period 1949—1953 showed 
that it is relatively simple to attain standardized conditions of counter operation with 
an efficiency unchanged from filling to filling. 

Ad. 2. Obviously neither the quality of the material and a. c. shielding nor the 
degree of the radioactive contamination of the materials of the counter depend on 
the technique employed. For this reason in discussing the techniques we do not 
consider the measured background levels published thus far by various laboratories. 
We assume that in principle the degree of contamination depends only on the size 


of the detector. 
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In order to determine the radioactivity of the samples two measurements are 
required: the number of background pulses in a unit time (T) and the number of 
pulses in a unit time of the test sample (L) + background. The standard error 6(L) 
of the end result for (Г) = (L + T) — (T) is 


6(L) = [9(L + Т)? + 6(T)2]*. (1) 


If we agree to regard as sufficiently accurate a measurement that is made over a speci- 
fied period of time (t) for which kó(L) <(L), (k > 1), then C = Lo", (Lo = Г/М), 
and M — the mass of the sample introduced into the detector for which a measure- 
ment is still possible with 6(Z,) < Lo 573) — can be defined as the sensitivity of the 
method. 

If the interval of time ¢ at our disposal was divided between measurements (L + T) 
and (T) in the proportion îr, т: tr then 
T--L n 


2 $ 
İT+L t іт Па) 


(Г) = 


For L> T and the application of the solid state carbon method using the S. W. C., 
ұт” tre 12; when a D. 5. W. C. is employed the measurement (T) and (T + Г) 
takes place simultanecusly, whence f, у= іт--і. When the gas sample method is 
used the background, in view of the fixed geometric conditions from measurement 
to measurement, may be determined once for a whole series of п measurements so 
that the actual time allotted for the measurement of the background for one measure- 
ment will be ¢/2n, the remaining time (2n—1)/n * t can be used for the measurement 
(T + L). In evaluating the value of Т and L for the various methods it is necessary 
to take into account the following: if the material and a. c. shielding are 100% effi- 
cent and the counting circuit is capable of eliminating double and delayed pulses,- 
then practically only the pulses originating in the radioactive impurities of the counter 
material and the material shielding take part in the production of the background. 
Generally speaking, this is therefore a surface effect, proportional to the surface enclos- 
ing the active volume of the counter. In this manner for the S. W. C. with an active 
length of lọ = sl and an efficiency of ns = 100% 


Ts = nfl? (ss + B[2) * Ty, (2) 


(where f is the ratio of the diameter to the length of the counter, ly, J the active length 
and the overall length of the counter respectively, Ty number of background pulses 
per cm? per minute). For gas counters operating at an efficiency of 7 


fg=c ВІ? CEDERA (2a) 


The number of 14C decays is proportional to the effective mass of radioactive carbon 
introduced into the active volume of the detector. 
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We therefore have 
ВИ, (3) 
and correspondingly, 


и. (3a) 


Lg = 7:90 Ру: Гос = прє: 


Е = 3.8% denotes the effective efficiency of an "infinitely" thick layer of carbon, 
о = 28 mg/cm? — the maximum range of 14C electrons, a= aq: w+p — the 
density of the carbon in the gas in g/cm? (до is the density of carbon under normal 
conditions in a gas in which the molecules contain one atom of carbon, w is the number 
of C atoms in a molecule of the gas sample, p is the pressure in atmospheres). Assum- 
ing l; = 15 = l and replacing in the formulae (1а) the appropriate values for L and T 
we find the ratio of the sensitivity of the counters as a function of their lengths: 


го б _ Los ,, б lc--ARs:DB)5 


EE c i © И 


3 
(Rs = s + 82; Re = (1—n7) no 1 20]; 
ume tege Tq CE pq nop SD ps (4) 


In comparing the sensitivities of the methods it is necessary to bear the following 
in mind: 

In the normally used counters f is of the order of 0.1 to 0.2, so that the active 
length of the gas counter is of the order of 0.6 to 0.8 of the overall length, in the case 
of the S. W. C. s — 0.3—0.2, so that the active length is 0.2 to 0.3 of the overall 
length, Е · о is about 1.1 mg/cm?; ж: o9: 5j: p is, for the average detector quality 
(п = 0.75, бұ = 0.535 mg/cm?; р = 1 atm.; w = 1), of the order of at least 0.4 mg/cm? 
and increases with the pressure and the number of C atoms in a molecule. 

In Fig. 1 the solid line represents the curve y = y(/). А counter filled with CO, 
under a pressure of 1 atm. was taken for comparison. The value T, was taken to-be 
5 x 1073 min”! cm. In the solid state carbon method we are compelled to make 
use of only the small amount of carbon at our disposal. The gas counter enables us 
to make better use of the sample. This can be of particular importance with larger 
detectors and in the case of small valuable samples which are destroyed during the 
measurement. Techniques employing gas samples make it possible to make measure- 
ments on carbon itroduced into the counter in the gaseous state under a reduced 
pressure p, of a quantity for which the condition Lp, > Горо is fulfilled. For most 
laboratories Ly is now of the order of 1—0.5 decays of !*C per minute per gram of 
carbon, i. e. of the order of several %, of the radiation of present day samples. This 
means that for present day samples the measurement can be made, generally speaking 
by using only 10% of the carbon normally introduced into the counter. In the case 
of the solid state carbon method the reduction of the material required for the measure- 

7 
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ment is very limited, if at all possible, in view of the standardization of the operating 
conditions. If one takes into account the more than 20-fold greater efficiency of gas 
counters it is seen that this method is far more superior from the standpoint of the 
economical use of the material. 

Fig. 1 (broken line) shows the mass necessary for making a measurement by 
the S. W. C. method as a function of the mass of carbon in the gas sample, with the 
same statistical accuracy, for both measurements obtained in the same period of 
time. 

These curves are derived on the assumption that both techniques make use 
of 100% of the carbon contained in the sample. In fact, however, in the elemental 
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Fig. 1. The change in the ratio of the sensitivity of the counters y — СвіСз as a function of the linear 

dimensions (solid line). The magnitude of the mass M s of carbon required for making the measurement 

in a given time with a given statistical accuracy as a function of the mass of carbon necessary in the 

case of the “gas” technique in the same time and with the same statistical accuracy (broken line). The 

relation between the linear dimensions of screen wall counters as a function of the dimensions of the 

counters in the *gas" technique capable of making the measurement with the same statistical accuracy 
in a given time. 


filling the counter. 
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carbon technique, for strict standardization, the thickness of the carbon must some- 
what exceed 28 mg/cm? so that the assumption that full use is made of the carbon in 
this method is not correct. Moreover, in view of the manner of preparing samples it 
is necessary to take into account the possibility that some carbon is lost during the 
chemical processes. 

In methods emploing carbon dioxide (CO, — proportional counter or CO, + 
+ CS, — С.М. counter) the question of loss of part of the sample material practi- 
cally does not occur. 

In the case of other methods, where gas samples (СН,, C,H,) are used, these 
losses are much less than in the case of the solid state carbon method. 

Ad. 3. With the solid state carbon method the counter must be evacuated each 
time and the layer of carbon replaced anew. This circumstance greatly complicates the 
work. Here it is necessary to standardize not only the process of filling the argon- 
-alcohol mixture but also, or even to a greater extent, the preparation of the carbon, 
the preparation of the paste, the application of the paste to the walls of the counter, 
the drying of the preparation, etc. 

Moreover, the carbon obtained by the reduction of CO, with magnesium is marked 
by strong adsorption,.which creates an additional danger of contamination with radio- 
active gases present in the air?. 

The difficulties cited above are substantially smaller or do not exist at all when 
gas samples are used. The entire process of burning and purifying the gas can take 
place ina vacuum apparatus, at least in the case of СО, being used to fill the counter, 
so that the contamination of the sample during preparation does not enter the picture 
at all. The change of samples takes place here simply by pumping out the gas. In the 
case of more complicated chemical processing (C,H, and СН.) there may be impuri- 
ties introduced together with reagents in contact with the atmospheric air, but the 
danger here is much less than in the preparation of paste from powdered carbon. 

Ad. 4. Control of the efficiency and operating conditions of the counter. Fluctu- 
ations in the background of the measuring counter. 

А. Variation in the degree of radioactive contamination of the material of the 
counter of a different level of intarsion of the impurities into the material used for 

In view of the fixed geometric operating conditions of the gas counters and the 
possibility of precise standardization of the physical-chemical processes for obtaining 
and purifying the gases in vacuum sealed apparatus, background fluctuations connec- 
ted with intarsion of the radioactive impurities into the counter walls or the sample 
tested cannot occur during the preparation and changing of samples. In the case of 
measurements of the activity of solid state carbon samples by the S. W. C. the pulses 
from radioactive impurities in the counter shield and the material used for filling are 


3 [t was found, for instance, that the degree of radioactive contamination of the carbon prepara- 
tions increased after the explosion of an atomic bomb in the Pacific (Blau, Deevey, & Gross 1953). 


218 W. Mościcki 


eliminated automatically, but only when the degree of contamination is the same in both 
samples (the inactive carbon and the sample tested). If the contamination is localised 
in any part of the apparatus or if the degree of contamination of the carbon prepara- 
tions is not the same, then the complete elimination of the background pulses does 
not seem possible and the results of the measurements will contain systematic errors. 

The D. S. W. C. have a certain advantage here over the S. W. C., but in this 
case also the possibility of misleading measurements due to local contamination of 
material must also be taken into account. 

B. Variation in the level of radioactive contamination in the laboratory and 
fluctuation of the cosmic radiation intensity during the measurements. 

In the S. W. C. a change in the background level induced by a rise in the radio- 
activity of the surroundings, by an increase in the intensity of cosmic radiation or 
possibly by a change in the efficiency of the a. c. shielding is noted during the counter 
operation and has an influence only on the statistical accuracy of the measurement. 

Distortion of the measurements because of background fluctuation does not here 
constitute any danger. 7 

When gas sample counters are used neither the size of the background nor апу 
fluctuations in it can be directly measured during the measurement. As a result, 
a very careful analysis must be made of the possible sources of background fluctu- 
ations in the measuring counter and of methods for controlling it. Changes in the 
number of background pulses in the measuring counter can come, when using gas 
samples, from three sources: a) changes in the intensity of у radiation of the labora- 
tory, b) changes in the efficiency of the measuring counter from measurement to 
measurement, c) changes in the effectiveness of the anti-coincidence shielding. In 
order to control the operation of the apparatus, pulses of the a. c. shielding, all the 
pulses of the measuring counter, and the anti-coincidence pulses of the counter 
shielding were recorded simultaneously during the runs with counters filled with CO, + 
+ CS. 

Let №, L and U designate the number of pulses registered in a unit time (cpm) 
for the anti-coincidence shielding, measuring counter and (L, N) anti-coincidences, 
respectively. To denote the origin of the pulse we use the following indices: k — pulses 
triggered by ionizing particles of the hard component of cosmic radiation; / — pulses 
due to y radiation from radioactive contamination in the laboratory; m — pulses 
practically not changing with time and due to contamination present in the walls of 
the counters and the shielding material; c — pulses from 14C electrons. All pulses 


in the shielding, measuring counter and the pulses И = (Г, N ) can be represented 
in the form of the sum: 


N= № + № + M, ay (Ny + М + Np), 
L= І, Же, + L, + a, (1, d Li Fl FL), 
LE P ЯБА. + xL, + ap (Lr Жеш tL. (A) 
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ау» ау, ау denote factors dependent on the operating voltages and the slope oft he 
counter characteristics, the time constants of the counting circuits or other parame- 
ters. The products ay N, a, Г, ay U denote the total number of pulses not directly 
connected with the passage of the ionizing particle through the detector or, in the case 
of the proportional counter, connected with the passage of a particle producing an 
initial ionization lying outside the range of the pulse amplitudes selected by the 
counting circuit. x L, denotes the pulses of the measuring counter produced by cosmic 
radiation and not accompanied by pulses from the anti-coincidence shielding. Such 
pulses may be connected with the inefficiency of the anti-coincidence shielding or 
possibly with the remainder of the non-ionizing component of cosmic radiation not 
absorbed by the shielding material. The means usually employed for eliminating y 
radiation (shielding material) and cosmic radiation (anti-coincidence shielding) and 
for ensuring the adequate stability of the work of counters and registering circuits 
make it possible to limit the number of pulses N,, Ly, ay М, aL and ар U to values 
that are small in comparison to the normal statistical fluctuations and fluctuations of 
cosmic radiation. The limiting of the factors for counting М, L and U effectively to 
three in number — cosmic radiation — (А), constant contamination (impurities in the 
material of the apparatus and the counter) — (m), and radiocarbon — (c), makes it 
possible to find convenient criteria for the proper functioning of the equipment from 
measurement to measurement and to detect any faulty operation. With sufficient 
stabilization of the voltage source the values of a can be maintained within limits that 
do not distort the measurement. For counters of the anti-coincidence shielding with 
a 295/100 plateau slope at the operating voltage and a counter voltage variation not 
greater than 2V, ay < 4X 1074, For М = 103 we have ay М < 0.4 cpm, which cor- 
responds to the standard error for a measurement lasting about 100 hours. In view of 
the small probability of uncontrolled fluactuation of the voltage in one direction during 
100 hours of operation the actual value ах, in an apparatus that is working well will 
not affect the statistical accuracy of the measurement of the number of pulses N. 
Similarly, a; L, as is shown by direct measurements, can be maintained within limits 
below the standard error for measurements lasting over several days. In using G. M. 
counters filled with a mixture of 250—400 mm Hg of CO, - 8 mm Hg of CS, operat- 
ing in the region of 50— 100 V above threshold the slope of the plateau does not exceed 
195/100. When L = 250 and AV = 10 V, then а, L < 0.25 cpm, which corresponds 
to a standard error for a measurement lasting about 3 days. If U «Є L, and if the statis- 
tical accuracy of the measurement U is to be better than L it is necessary, when record- 
ing U, to employ special means to ensure that a, U < (U) < a; Г. For this reason 
the elimination of counts due to spurious pulses in recording U must be done with 
particular care. Adequate stabilization of the voltage supply and elimination of spu- 
rious pulses from the counts allow the factors to be made sufficiently small so that 
the corresponding products do not place a limit on the statistical accuracy. 

The separation of the N, Г, and U pulses into components depends on the quality 


of the material and anti-coincidence shielding, and. on the manner of location of the 
j 
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counters of the a. c. shielding and the measuring counter. With a particularly good 
shielding U can be of the order of 1.5 to 2%, L (Kulp & Tryon 1952). The characteris- 
tics published for various types of equipment show that shielding of 150 to 200 gm/cm? 
ensures, practically speaking, complete elimination of the soft component of cosmic 
radiation and y rays of radioactive contamination in the laboratory. During the work 
with G. M. counters filled with CO, + CS, with a shielding of 10 cm of lead it was 
found that the efficiency of the elimination of y rays from radium in equilibrium 
with decay products was not worse than 99.7 = 0.2%. In one brass counter (— 2.5 2) 
filled with a mixture of CO, 4- CS,, used in this work there were noted in an unshiel- 
ded counter about 600 cpm due to radioactive impurities in the laboratory. Inside the 
shield, therefore, 1.8 to 1.2 cpm were to be expected from this source. In this counter, 
which was placed inside a shielding, it was found that Ly ~ 250 cpm and U, ~ 80 cpm. 
When the brass counter was replaced by a steel counter of almost the identical dimen- 
sions and the material shielding separating the counter from the lead shielding incre- 
ased by 6 mm of Hg the number of pulses dropped to Ly ~ 200 cpm and U, ~ 35 cpm. 
In these apparatuses the probable values of the L component were L, ~ 170, 
І, ~2, La ~ 78, in the case of the brass counter and Г, ~ 170, І, — 2, and І, ~ 33 
for the steel counter. 

Іп опе of the 'cleanest' apparatuses (Kulp & Tryon 1952) Г, of the unshielded 
counter was about 225; with a shield of over 20 cm of Fe an L, of the order of 0.2cpm 
is to be expected inside the shielded counter. The probable pulse distribution in this 
apparatus was, for the counter of ordinary brass 


о 


and for a counter of a specially selected material surrounded by distilled mercury and 
with. an increased anti-coincidence shielding 


PR o. eu d W ml LO apd. ро ау: 


From the examples cited it can be seen that in the case of the normally used 
shielding (150—200 gm/cm?) changes in the background level connected with 
fluctuations in the intensity of the laboratory radiation of a composition near that 
of Ra radiation in equilibrium with the decay products, will not limit the accuracy 
of the measurement if the amplitude of these fluctuations does not exceed 5—10% 
of the mean intensity. Probable changes in the background level connected with 
this source can, in principle, be assessed from the control measurements of the 
laboratory radiation intensity and from a knowledge of the efficiency of the 
material shielding. In work with СО, + CS, counters (1950—1953) no fluctu- 
ations, of long duration, in the laboratory y radiation exceeding 5% of the mean 
intensity were observed: From the examples cited above it can be seen that the Г, 
component constitutes the main part of the L pulses and Г, of the U pulses. In any 
event, even for apparatuses of weaker shielding L, > Lp > Г. If the material shielding 
is good enough to ensure the elimination of diae feni the radioactive contami- 
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nation in the laboratory so that the uncontrolled changes in the number of background 
pulses connected with changes in the level of the laboratory radioactive contamination 
AL, < 6(L,,) and if at the same time x ~ 0, then in view of the fixed geometry of the 
gas counters, the factors responsible for the pulses are reduced. Instead of (A) we can 
write: 


N= № падаў + М), 
Li tis “ei a dc 
AU aoe 


(N,, + №) and (Z, + Г.) remain constant (within the limits of the statistical fluctu- 
ations) from measurement to measurement. The changes in N and L—U observed 
from measurement to measurement can come (with an apparatus that works well) 
only from changes in the intensity of the hard component of cosmic radiation. If, by 
chance, for two different measurements М, = N, and (L—U), = (L— U), within the 
limits of the probable statistical fluctuations N and L, then there is no reason to doubt 
that the operating conditions for both apparatuses were identical. 

If, however, AN = N,—N, = AN, > 6(AN) and A(L—U) = AL, > 6(AL,), 
then the ratio k = (4М-- A(L— U)/ A(L— U) should be constant, within the limits of 
the statistical fluctuations. 

The accurate determination of the statistical distribution of k does not seem to be 
possible. An approximate distribution can be obtained by introducing certain simpli- 
fying assumptions. ; 

Among the anti-coincidence shielding pulses produced by the hard component of 


cosmic radiation — М, a part, №, = М, — L, is statistically independent of the Z, pulses. 


C, = №, Г, is an equipment constant; it does not change from measurement to meas- 
urement and is independent of the momentary values L, and N,. In reality, for a speci- 


fic measurement N,/L, = С >< Со. Let us first consider the statistical distribution of 
the ratio C. For the sake of simplicity let us assume: 1. During each measurement the 
intensity of the cosmic radiation remains constant. 2. The statistical distribution of 


the number of pulses N, and Г, can be approximated by a Gaussian curve. 
In accordance with the assumptions the probability of counting a number of 


pulses N, contained in the interval y, y + dy and the number Г, in the interval x, 
х + dx will be 


С (Муз у) dy = (2л oy) ^ exp [— (y — №32 of] dy 
and | | 2 
С(Г,,х)4х = (2л ор) ^ exp [— (x — 1,)92 02] ах, 
and since the Gaussian curve is here an approximation of the Poisson distribution 


then 02/02 = ҚАТА = Со. The probability that the ratio y/x will be contained within 
the limits C, C + dC will be 
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malin (s — Ly) (Cx — Мы 
Р(Со C) dC = (дл огом)? ШЕЛ e Aor бе? eye | dC 
0  —oo 


"Cod Cu) (дл) (Cot Сб» 2% (Cot C?) 


x 
(ô = Г.с) and © (x) = 2/л% thos dt). 
б 

With the statistical accuracy of the measurements of the activity of natural carbon д 
is of the order of 10?— 103, so that the first term of the expression becomes small, 
Ф(х) - 1, and the entire second term differs distinctly from zero only for values of С 
close to Cy. This makes it possible to replace the distribution P(C,, С) by the Gaus- 
sian curve distribution 


С (Cy, C) dC = (2л а) "Я exp (— (C, — 042 o?) dC 


where 
e, = [Со (C, + 1)%/6|. 


Since C, is usually of the order of several units o, < 10-2 and the chances of obtaining, 
during the measurement, a value of C differing by more than 1% from C, is usually 
negligibly small. 

It therefore seems reasonable to expect that for each pair of measurements for 
which the expression is determined there must hold the inequation k— Ak < C, < 
<k + Ak, where Ak denotes the total inaccuracy in the evaluation of k and is con- 
nected with the standard errors of the measurements of AN and A(L— U) and other 
determinable factors or errors. 

In order to get an insight into the distribution for the probability of obtaining, 


for a given pair of measurements, values for the ratio AN,|AL, contained within the 
limits k, k + dk it seems justifiable to assume within the framework of the previous 
assumptions a distribution P(C,, k) where we introduce n = AN, in place of N, 
l = AL, in place of L,, 6, = в, o, = [02 + 02 |* and o, = loi, + ox] in place of 
д, от, and oy resp. The shapes of the distributions for С, = 5 and various values 
of the parameter д, are shown in Fig. 2. 

Owing to the lack of suitable data in papers published thus far on apparatuses 
for which the assumed conditions are in certain fulfilled we shall employ for the sake of 
illustration data obtained from the author's experimental apparatus during the period 
1948 to 1950 in Poznań,. 

The counts obtained for various samples are given in Table I (Moscicki 1953). 

The table contains the following data: col. 1 — sample symbol; col. 2 — date 
surement; col. 3 — duration of measurement (in min.); col. 4 — counts per min. for 
of meathe a. c. shielding; col. 5 — total number of pulses in the measuring counter; 
col. 6 — number of anti-coincidence pulses U; col. 7-— number of coincidences pulses 


) 
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(L— U); 8, 9, 10, 11 — standard error N, Г, U, (L—U); 12 and 13 — corrections 


for instability of the operating of the counters; 14and 15 — maximum error in the 
measurement of N and (L—U). 
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Fig. 2. The distribution function P for various values of the parameter 6 and the constant Со = 5. 


Table I 
Good 
B Good 
c 3.60|0.43| Good 
D 2.99|0.19|4.09|0.52| ~85% effi- 
{ сїепсу 
Е 2.75) 0.23|4.85|0.93| Hole in shiel- 
ding (Duea.c. 
counter out) 
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Table I gives the values of k for the various pairs of measurements and the limits 
of the values of k resulting from the total possible error in measurement (Table I, 
columns 14, 15). In parentheses are given the values of д for the various measurements. 


Table I 
he pec | DU p 
| >= 
4.83 2.10 —24 — 80 
А 21.40 72.90 0.9 5292 
(6 = 19) | (6 = 4.2) К 
6.82 — 18.6 0.6 
B -=2.80 | F68 +13 
(б = 6) 
ПЕ 49 TUE 
C +15 Е 5.5 
Са 
D +12 


For measurements A, В and C, which were found to be good, the following limits 
for k are obtained by taking into account only the standard error: 


5.2 > k (AB) > 4.5; 29 > ЩАС) >17; 7.7 > (BC) > 6.2. 


As may be seen from the listed data not all of the assumptions necessary for main- 
taining the statistical constancy of the ratio k were maintained in these measurements. 
The a. c. shielding of this apparatus, it is true, did indeed consist of 21 counters with 
plateau thresholds within the limits of about 150 V., the slope ofthe plateaus of the va- 
rious counters varied within limits of 2—5%/100 V so that the slope of the plateau 
of entire anti-coincidence shielding was, at the operating point, about 595/100 У. 
The voltage applied to the separating bridge was controlled to an accuracy of 10 V. 
Under these conditions ay was of the order of 5 x 10-3. The L pulses were registered 
directly on a scalar with a resolving time of 10-5 sec., so that there was no elimination 
of possible afterdischarge and multiple pulses. The observed slope of the plateau of 
the L pulses was 1—2%/100 V in the operating region of the counter. The maximum 
difference in the counter's high tension readings in two different measurements did 
not exceed 20 V. Under these conditions ау could attain to 2 x 10-3. Taking into 
account the maximum corrections (Table I) connected with possible differences of 
voltages of the operating points of the counters we obtain: 


k(AB) = 4.83 1.40 
ҚАС) = 2.10 2.90 
(ВС) = 6.82 = 2.80 


These results do not allow us to draw negative conclusions аз to the operation of the 
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equipment in any of the above three measurements. For comparison, in Table I and 
П lists two measurements made with incorrectly functioning equipment. Measure- 
ment D refers to a counter filled with CO, -- CS, with an admixture of air. The 
efficiency of this counter was estimated to be at about 85%. In measurement E one 
of the 14 side counters of the anti-coincidence shielding was disconnected. 

The statistical stability criterion of the quotient AN—A(L—U)/A(L—U) is 
particularly sensitive to the difference in efficiency of measuring counter in two 
comparative measurements. Denoting the relative efficiency of one the counters with 

(14-À)e 
e+ (y—1) 
where ғ denotes the relative increase in the intensity of the cosmic radiation registered 
by the anti-coincidence shielding e = AN,/N,. Usually e does not exceed 0.1, and 
hence even very small differences of the measuring counter efficiency can result 
an enormous increase of y, and when y < 1--є, y becomes negative. 

А d. 5. As has been stressed repeatedly (Anderson, Levi) the electronic equipment 
employed with G. M. counters is much simpler and easier ot use than in the case of 


respect to the other by y we have for the registered quotient the value y — 


1, 


proportional counters containing gas under a pressure of many atmospheres and 
requiring high (of the order of several kV), well-stabilized voltages. The apparatus is 
complicated by the necessity of using linear amplifiers and pulse amplitude selectors. 
The power supply and registering equipment in apparatus using the G. M. counters 
filled with a mixture of CO, + CS, are much less complicated than in other equipment 
employing the gas sample method, but in view of the high voltages required and the 
necessity of triple registration are more complicated than in the case of'the solid 
carbon technique. 

Ad. 6. It is evident from the above considerations concerning the sensitivity of 
equipment for the measurement of the activity of 14C that the sensitivity is a function 
of the increasing geometric dimensions of the counter. Ап increase in the geometric 
dimensions of the counter entails the necessity of increasing the material shielding and 
increasing the overall costs of the apparatus. As has been mentioned in section 1, 
the gas sample method has an advantage in sensitivity, particularly when large detec- 
tors are used. Fig. 1 shows the direct relation between the linear dimensions of the 
. counters of the S. W. С. type and the gas counters filled with a gas having ап atom 
of C per molecule at a pressure of 1 atmosphere and capable of making measurements 
with the same statistical accuracy in the same time. 

Moreover, considerations of cost and size of apparatus employing counters one 
metre in length call for rational design of material shielding and anti-coincidence 
shielding. Material shielding in the form of a *barrel* and an anti-coincidence ring 
(Jóźwiak and Mościcki 1953) with the addition of a “гооҒ“ ої anti-coińcidence counters 
on top seem to be the most economic solution as regards the quantity of material 
used and space required. « 

In equipment of the S. W. C. type this solution, in wiew of the continual necessity 
of removing the counters, does not seem possible. 


i 
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КРАТКОЕ СОДЕРЖАНИЕ 


В. Мостицкий, Опріїменений счётчика Г.М. наполненного смесью СО, + CS, для йзмеренйя 
активности натурального угля. 

Часть I. Сравнивается технику измерения радиоактивности натурального 
угля а) пользующейся углем в сталом виде и счётчиками S.W.C., б) пользующейся 
газовыми пробами. Техника газовых проб обладает превосходством по чувстви- 
тельности израсходования материала, вымене проб, точности очистки. Трудно- 
стью техники в случае применения пропорционального счётчика является слож- 
ная аппаратура, а особенно отсутствие контроля уровня шумов во время изме- 
рений. Обсуждается статистический метод контролирования шумов использую- 
щих флюктуацию интенсивности космического излучения в счётчике и перего- 
родке. 


і 37) 8 ў ^ P M wa. аяғы ви 


о 32 жд раба іо 
(iari naa pe. A ET де де age E, 
a te P 


C—— P 


Vol. XVII (1958) - ACTA PHYSICA POLONICA Faces 


ON THE USE OF G. M. COUNTERS FILLED WITH A MIXTURE OF 
СО, + CS, FOR THE MEASUREMENT OF THE ACTIVITY OF 
NATURAL CARBON 


Ву W. Mościcki 


Part II 


EXPERIMENTS WITH С. M. COUNTERS FILLED WITH СО,-- CS, USED FOR DETECTION OF “С 
IN NATURAL CARBON 


(Received march 13, 1958) 


A number of experiments were conducted with С. M. counters filled with CO, + 
+ CS,. The counter has a long plateau with a small slope and low sensitivity to contami- 
nation. The change in the number of the counter-shielding anticoincidences as a function 
of the delay time of the shielding signals is interpreted as resulting from pulse delay in 
counters of this type. According to this interpretation about 65% of the pulses are subject 
to a mean delay of about 0.8 m sec. The influence of the delay effect on the efficiency 
and suitability of the counter for chronometric purposes are also discussed. Reproducible 
results are obtained with a sufficiently long time delay of the a. c. shielding pulses 
(~ 5 пі вес.) over a broad interval of СО, pressure and counting efficiency of about 
90% for counters with a volume of 2—3 litres. 


1. Introduciion 


It is known that pure CO, cannot be used to fill G. M. counters because of the 
considerable probability (10-4) of negative ions being formed upon ion-cathode соШ- 
sions. Müller (1947) showed that the addition of several % of CS, produces a mixture 
suitable for G. M. counters. Such counters require external quenching. Müller studied 
the properties of counters of this type for small counters and partial pressures of CO; 
of the order of several tens of cm of Hg with the addition of several % of CS. 
Miiller‘s counters exhibited very good and long characteristics and were employed 
in work on determining the half life of 14С. In order to obtain reproducible results 
сіп comparing the results of the activity measurements of samples of different concen- 
trations of МС it was necessary to introduce a “linearity correction factor“ by assuming 
values from 1 to 1.295. Depending on the shape and filling mixture the efficiency for 
counters of various shapes was estimated at 95— 10094, and for counters of a given 
shape the output was entirely independent of the CO, pressure. The slope of the plateau 
of these counters did not exceed 395/100 V. The work of Hawkings et al. (1948) and 
Mann and Perkinson (1949) seemed to indicate an efficiency of 100%, for counters of 


+ 
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this type. A direct comparison of the values of T4 for 14С obtained by various methods 
including the CO, + CS, counter showed that for unexplained reasons the efficiency 
of this counter is about 85%, (Miiller et al. 1950). 

From the point of view of application of counters of this type for chronometric 
purposes a very important circumstance is the above-mentioned fact that the efficiency 
of counters of a given shape is constant and independent 
of the CO, pressure. The aim of this work was to study 
the conditions for reproductobility of the results in the 
measurement of !4C radiation from natural carbon. 


2. Apparatus 


a) Counters. А. с. shielding. 

In the first experimental apparatus! a counter de- 
signed by Miesowicz and Jurkewicz (1947) was used аз 
the measuring counter (Fig. la, b). For the anticoincidence 
shielding were used counters of similar construction with 
a length of / =300 mm, anode length / = 270 mm, dia- 
meter Ф = 30 mm for lateral shielding, and / —155mm, 
1, = 130 mm for the top shielding. Anodes consisted of 
Kanthal wire Ф = 0.1 mm.The counter was filled with a 
mixture of 90 mm Hg Ar + 10 mm Hg ethyl alcohol. The 
plateau thresholds of the a. c. shielding counters were 
within the limits of 150 V.The length of the plateau was 


Fig. la Fig. 1b 


Fig. la. Design of the measuring counter, а — brass tube and lid, b — brass-glass couplings, c — steel 
springs 2 0.25 mm, d — glass insulators, e — anode wire. 
Fig. 1b. Glas insulator couplings with brass. a — anode wire, b — layer of silver, c — electrolytic copper. 
і d — tin, e — shellac. 


150 to 250 V and the slopes of the plateaus were within the limits of 2 to 595/100 V. 
The a. c. shielding consisted of 21 counters of which 14 were lateral counters (2 rings 
of 7 counters). In later work (Torun, 1952—53, Gdańsk 1955—56) counters of the 


! Built in connection with the author's doctor's thesis in the Department of Mathematics, Physics, 
and Chemistry of the University of Poznań, 1951, “On the Measurement of the Concentration of the 14C 
Isotope in Organic Carbon.” 
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А- 


design shown іп Fig. 2 а, b, с were employed. The upper a. c. shielding in this case 
consists of 14 counters of / = 300 mm, J = 270 mm filled with an argon-alcohol 
mixture and placed directly above the measuring counter. The characteristics of the 


400 
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Fig. 2a. Design of the large measuring counter. a — rubber or igelite seal, b — leads from the ring 
electrodes of the a. c. shielding, c — screw cap sealing the opening of the glass insulator, d — tin, 
e — shellac. 

Fig. 2b. Leads from the ring anodes of the a. c. shielding a — seal, b — araldite, c — shellac. 
Fig. 2c. Manner of fastening the ends of the electrodes of the ring shielding. 


. counters used are given in Table I. The set-up of the counters and shielding in both 
cases are illustrated in Figs. 3 a, b and 4 a, b. 
/ 
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Table I 
Total Active Cathode Anode Background 
Counter | volume volume Material © Length Material Ø (cpm) 
(стз) | (em?) (ст) (cm) (cm) 
СТУ er 

N 23 432 385 Silver 4.5 24.2 Kanthal 0.01 10.7 

М 37 439 360 Brass 4.5 22.1 Tungsten 0.01 11.2 
M2 2940 1810 Brass 9.8 24.0 Kanthal 0.015 80 
MA 3145 1940 Steel 10.2 23.8 Kanthal 0.015 35 


Fig. 3. Location of the anticoincidence counters (4) of the lateral shielding, a) in the first work, b) in the 
second work, a — anodes, k — cathodes of the ring shielding, Hg — 6—10 mm layer of mercury. 


b) The materiał y — shielding was comprised of a “lead barrel“ 40 cm high 
with 10 cm walls (Figs. 4a and 46). In the later equipment there was an additional 
mercury filter first placed between the lead and the anti-coincidence ring and after- 
wards between the measuring counter and the a. c. ring (Figs. 2, 4b). 

Direct measurements of the efficiency of the material shielding by comparing the 
increase in the number of pulses of the unschielded counter and inside the shielding 
irradiated by a radium preparation with an unchanged geometric configuration, showed 
that the efficiency of the shielding for the absorption of y rays from Ra in equili- 
brium with the decay products amounts to 99.7 Е 0.2% (within the limits of one 
standard error of measurement). 

c) Electronic Apparatus 

The block diagram of the electronic equipment is shown in F ig. 5. The h. t. 
supply allowed continuous adjustment of the voltage within the range of 0—12,000 V. 
For technical considerations the maximum voltage that could be used for the measuring 
counter was not over 5 kV. The voltage was stabilized by means of an A. C. stabilizer 
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Fig. 4. Cross section of the material shielding, counter shielding and measuring counter a) in the first 
work, b) in later work: a; — lateral ring shielding, a — counters of the upper shielding, c — measuring 
counter, Hg—6—10 mm of mercury, Pb—10 cm lead “barrel”, Fe—table frame 10 cm wide, З cm high. 


Fig. 5. Block diagram of the electronic equipment: H. V. — high voltage supply, continuously adjustable 
over the range 600—12.000 У” V. S. — high voltage stabilizer” N. P. — quenching circuit” К. О. — 
cathode oscilloscope; R — registering apparatus (Scaler 1:64) for pulses N; R, registering apparatus 
(inonostable univibrator with a delay time of 107? sec) of pulses U; Rg — recorder (scaler 1:8) of pulses L; 
P — stabilized supply of the shielding counters; S. C. — anodes of the counters of the a. c. shielding ; 
KF — cathode followers; 4C — anticoincidence system with an adjustable multivibrator time delay 
| for the a. c. channel. 
4 i 
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(Grzesiak 1954) and an ,,s" type stabilizer (Hount & Hickmann 1939) with a cathode 
resistor (Mościcki 1951) having a stabilization factor of about 200. Apart from the vol- 
tage stabilization, there was a constant control of the counter voltage by measuring 
the pulse amplitude with the use of a cathode oscilloscope. 

То quench the measuring counter a modified Neher & Pickering system (1938) 
was employed — Fig. 6 — and in the later work (Тогай 1954) a special power supply 
and quenching circuit (Gorgolewski 1953) was used. 

In view of the different threshold voltages of the a. c. shielding counters, the 
voltage from the shielding counter supply was fed through a divider bridge. 

In this work the anticoincidence circuit was of the du Toit type (1947) in which 
the widths of the multivibrator gates can be smoothly adjusted over the entire range: 
counter pulses from 50—300 u sec and anticoincidence shielding 1—10 m sec. The 
a. c. shielding pulses (N) were counted by a 1:64 scaler, pulses of the measuring 
counter (1) by a 1:8 scaler, pulses (U) of the (Z, М) anticoincidences by a circuit with 


Fig. 6. Quenching circuit. 


a constant blocking time of 107? sec. In work with small counters (Poznań 1950—51) 
the pulses U were recorded directly by a telephone counter, in later work with big 
counters by means of a 1:4 scaler. 

d) Combustion of Samples, Purification and СО, + CS, Dosing. 

The block diagram of the apparatus is shown in Fig. 7. In the first work with 
small counters organic samples were burned in a vacuum «sealed apparatus by means of 
a van Slyke-Folch mixture (1940). The burning of such samples as wood, peat, char- 
coal, horn, amber, etc. was done cold. In the combustion of bituminous coal and 
anthracite it was necessary, for starting the carbon oxidation, to heat the mixture quite 
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strongly (80—100°C). The carbonates were decomposed by means of 10% НСІ. In 
later work the samples were burned in an atmosphere of oxygen in a combustion 
tube. In these cases the СО, obtained was combined by means of freshly prepared 
NH,OH, then changed into СаСО., by means of CaCl. To obtain CO, for the counter 
the process was continned in a vacuum apparatus. CaCO, was treated with 10% HCl 
and later (Gdansk) with an aqueous solution of citric acid. 


The V.S. F. mixture and a fine sample of the organic substance (of a mass of 
2—5 gm) rinsed with 10% HCl and then with distilled water and dried, were intro- 


Fig. 7. Block circuit of the apparatus for the combustion, purification and dosing of CO, + (52, 

A — V. S. F. mixture; S — container for burning samples, 46 — filter of silver wire T, — Т, freezers, 

CaO — quartz tube with calcium oxide R, — Ёз — containers, CS, — container for carbon bisulphate, 

P — mercury pump, С — counter, M — М, — manometers, V — vacuum, Ку — Ку; — vacuum valves, 
Zı — Z, — mercury vacuum valves. 


duced into the container S. Then all of the valves were opened and the apparatus was 
pumped. down to a pressure of 10-1 mm Hg. Valve K, was opened and the V. 5. Е. 
mixture was allowed to flow through the condenser into the container S. About 250 cm? 
sufficed to burn the samples for the small counter; in work with the large counter the 
combustion process took place in two stages. In this case about 1.5 | оЁ mixture was 
used. The process of combustion, depending on the type of sample, lasted from 
30—90 т. The CO, produced as the result of the combustion was passed through a tam- 
pon of thin silver wire heated by a small flame. Then it was subjected to repeated distilla- 
tion at a temperature of about —114°C. АП of the gas was collected in the trap Ту at the 
temperature of liquid air, the valve Ка was closed and after careful emptying of the 
E i 
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apparatus, was distilled through valve К, at the freezing point of alcohol (—114.15°С). 
into the freezer cooled by liquid air. In a similar way the CO; was distilled from T; 
into T, and from T, into Ту. The thermal insulation of the freezers was such that the 
rate of evaporation of the carbon dioxide was almost constant during the evaporation 
process and amounted to 20—30 cm? a minute (under standard conditions). The 
first batch of the gas evaporating from the first freezer was discarded (50—100 cm? 
of CO, under standard conditions), the remainder being collected in the next freezer. 
The CO, evaporation process was never completed. When only slight amounts of 
CO, remained in the freezer the freezer was cut off and the process of evaporation was 
stopped. The CO, losses in this method of purification amounted to 19—20%. 

The repeated distillation at a temperature of about —114°C makes it possible 
to reduce to a minimum the effect of the mechanical carrying of water or dilute acids. 
In some cases the total quantity of CO, impurities that got into the first freezer amoun- 
ted to a volume of 0.5—1 cm?, and hence in the evaporation traces of impurities 
were carried off mechanically by the CO, stream and were carried into the next freezer. 
The sublimation temperature of CO, (—114°C, the freezing point of alcohol) was 
used in order to reduce to a minimum the pressure of foreign substances. At the 
temperature commonly used for evaporating СО, under normal pressure (—80°С) 
these substances still have a considerable saturated vapour pressure. For instance, 
$05, which is likely to be present, especially in the burning of anthracite, has a satu- 
rated vapour pressure of about 8 mm Hg at —80°C. The ratio p,, /р,, at this tempera- 
ture is 95, whereas at —114?C it is about? 450. 

The saturated vapour pressure of CO, is, at —114^C, still high enough 
(— 25 mm Hg) so that given a large sublimation surface, the distillation does not take 
very much time. With this method of purification one can be certain that at the 
pressure of CO, used (about 0.5 atm) the partial pressure of SO, in any case never 
exceeded 1 mm of Hg. Apart from this method of purification there was also employed 
a modification in which there was a quartz tube with pure CaO between freezers T, 
and Т,. The CO, from freezer T, maintained at about —114°C evaporated into the 
quartz tube containing CaO where, at a temperature of about 600 to. 650?C, it was 
combined into СаСО. Once all the gas had been combined the valve was closed, 
the temperature of the oven was brought down to about 400?C and the pressure was 
pumped down to 10-2 to 10-3 mm Hg. After 11], to 2 hours of pumping the valve 
was closed, the oven temperature raised to 800—850?C, and the resulting gas was 
collected in the freezer from which the CO, was either slowly evaporated through the 
freezer (kept at about —80°С) or distilled in conditions as above to T, and from there 
into the container. The last two freezers, moreover, enabled the convenient manipu- 
lation of quantities of CO, for producing the required mixture of CO, + CS,. The 


? Unfortunately, is has not been possible to find the pressure of SO, at —114? C in any available 
literature. In evaluating the ratio р, [р;,, ghe value p(SO;) was evaluated by extrapolation of the experi- 
mental curve р = p(t). he 
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СО, was introduced into the container by the intermediary of ап auxiliary container 
(Fig. 6). Since, before the CO, is let into the container R,, this container already 
held a *standard* portion of CS, it was necessary to introduce the CO, in such a way 
so that the pressure іп R, was always larger than in R,. The same chemically pure 
С5,, vacuum sealed in the (CS,)container was used throughout this work. The desirable 
quantity of CS, that would give a pressure of 8 mm Hg in the measuring counter was 
obtained by cooling the CS, container and the auxiliary container R; of appropriate 
volume to 0°C, the apparatus pressure being reduced to 10-3— 10-* mm Hg. Container 
К. was filled with saturated vapour of CS, at 0°C, the flow of CS, was stopped, and 
CS, was let into Ry. 

The СО, pressure in the container was so chosen that after 2/3 of the gas had been 
pumped into the counter by means of a mercury pump (P), the desired pressure of 
the mixture was obtained. The magnitude of tlie final pressures in the counter and the 
container was measured on manometers M, and M,. 

The entire operation of combustion of the preparation and the purification ofthe 
CO, took about 12 to 18 hours and was usually conducted in two stages: 1) combustion 
and combination of the gas in the furnace into Са СО. ; 2) production of CO, by the dec- 
omposition of CaCO, and the introduction into R,. No systematic differences related to 
the manner of combustion or method of purification were observed in the behaviour of 


the counter filled with the mixture of CO, + CS,. 


3. Some Properties of the СО, + CS, Counter 


a) Thresholds and Plateau Length. 

In view of the purpose of this work — to find the conditions for the maximum 
standarization of the counter operation — the material contained in previous papers 
(Müller 1947, Brown & Müller 1947, Mann & Parkinson 1949) and insufficient techni- 
cal equipment in the laboratory, no quantitative studies were made into the relation 
between the composition of the mixture employed and the threshold voltages, the 
slope and lenght of the plateau. Table II contains data concerning some measurements. 
Fig. 8 shows the typical counter plateau. 

b) Delay of the CO, + CS, counter pulses. 

It was observed that within the region P of the ges (Fig. 8), with the 
use of a Neher-Pickering quenching circuit, the number of pulses registered depended 
on the sensitivity and resolving time of the registering equipment. The typical shape 
of the majority of the pulses observed on the oscilloscope screen in this region is 
shown in Fig. 9. 

It was found that the values of the paramaters of the N. P. circuit do not have 
any influence on the time intervals between the various maxima (of the order of several 
tens of miliseconds). Nor was there any distinct influence of the СО, molecular pres- 
sure in the counter on the time intervals at which the maxima appear in the range of 
pressures used in this „work. 
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Table II 
= 
Total Threshold Length of Optimum slope of plateau 
Counter Type Pressure (start of plateau**) in % —100 V 
(mm Hg) plateau) kV Sr Sy? 
A | ES m" 

N 23 G. M. 250 2.9 200 «296 «196 
N 23 G. M. 278 3.35 180 2% <1% 
N 23 G. M. 289 3.36 200 <2% <1% 
М 37 С. М. 368 3.9 220 <2% ~0.5% 
N 37 G. M. 375 4.0 210 205 ^ 0.596 
M2 G. M. 400 4.65 260855) «208 —0.5% 
М4 С. М. 390 4.5 300 <<3% ~1% 
М4 С. М. 500 4.95 100553) E4900 ~1.5% 
MA Prop**) 320 3.55) 700 —5% 
MA x 350 3.5+) 700 --5% 
М4 ж 380 3:01) ы OD ~3% 
M 4 5 400 засц) 1000 ~3% 
Ма yi 500 4.27) 800***) ~1% 
М4 " 650 4.6*) 500***) ~1% 


*) With zy = 8.33 x 10-5 min. In the area A B (Fig. 8). 
жж) From the.threshold to the voltage at which L(V) = 1.2 LIV opt. 
***) In view of the limited supply voltage. à 
+) With discrimination imp. < 50 mV the voltage for which М = 0.9 N opt. 
++) CO, for filling counters operating within the area of proportionality was obtained by chemical 
decomposition of pure calcium carbonate. The CO, obtained, after combination into СаСО; chemical 
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Fig. 8. Typical characteristic of а С. M. counter filled with a mixture of CO, + С5,. The characteristic 
- shown is for a counter, filled with about 250 mm Hg СО, and 8 mm Hg С). P. A. denotes the beginning 
of the Geiger area. The shaded area is the opitimum operating region for the counter. 
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Outside this region multiple 
pulses were no longer observed. The 
statistics of time intervals between 
the pulses, for a counter with various 
fillings also did not show more than 
1:10? multiple pulses, within the 
limits of statistical error. 

When thelength of the signal 
of thea. c. shielding pulses was chan- 
ged a monotonic drop was observed 
in the number of pulses U of the 
(L, N) anticoincidences. Table III 
gives the data obtained for two 
different counters. 

The first eight positions refer to 
the large counter (Table I) filled with 
a mixture of CO, + CS, under a pres- 
sure of 400 mm Hg. Positions 9 to 16 
refer to the small counter /V 37 filled 
with a mixture of 376 mm Hg 
CO, --8 Hg CS,. The vertilac co- 
lumns give: 1. time of duration (Ty) 
of the pulse of the a. c. shielding 
closing the apparatus; 2. the registe- 


3 msec 


"AGS 


red number of pulses №; 3. the nu- 
mber of pulses of the a. c. shielding 
counters № = М (1--Лт,)-1; 4. the 
number of pulses of the measuring 
counter [;5. the recorded number 
of (L, N) anticoincidence pulses; 6. 
the total number of anticoincidence 
pulses of the measuring counter 
Л = U (1—Nry—Ut;);" 7. the 
ratio AU,/L in percentages, where 
AU, = U,(t)—U, (co) and 0, (со) < I о с 
denotes the mean value obtained » >, " 
from very large statistical material, 
where Ty = 8.33 х 10-5 m. The value 
of U,(cc) was 97.92 for the large counter and 13.73 for the small one. 

Fig. 10 shows the experimental points according to table V. The solid line repre- 
sents the function Z = 0.645 exp (—0.72 x 105/min) that best agrees with the 
experimental.data. As can be seen from the curve the experimental points for both 


Fig. 9. The shape of the pulses in the transitional region. 
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Table III 


U.(t)— (со), 
PIE Я ЗЛУ 


U, cpm 

147.7 - dE 

126.0 11.2 

110.0 4.0 
99.2 0.6 
93.5 —1.8 
98.0 0.0 
95.7 —0.9 
97:92” 0 
29.4 18,9 
18.0 4.95 
15.9 2.57 
15.6 ^ 2.58 
14.3 0.78 
14.6 1.09 
13.5 --0.24 
13.73 0 
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Fig. 10. The distribution of the delays of the pulses from the G. M. counter filled with a CO, 4- CS, 
mixture. On the ordinate axis is the % of delayed pulses; on the abcissa, the delay time. | 


counters fit well to the curve of Fig. 10. In view of the low accuracy (about 10%) 
in determining the duration of the anticoincidence shielding pulses and the delay 
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with which the measuring counter pulses appear — about 0.5 x 10-5 min in both 
cases — and the statistical character of the counts, the parameters of this curve are 
found with an accuracy of no better than about 20%. 

The function U = (т) seems to indicate that a considerable part of the pulses 
from the measuring counter are subject to a delay in relation to the moment at which 
the ionizing particle passes. For t — 0; Z(0) — 0.654, from which it would appear 
that about 65% of all pulses L are delayed with a mean delay of about 1.4 x; 10-5 min— 
= 0.84 x 10-3 sec. The mechanism of the pulse delay was not studied. If the time of 
flight of the ions is т, and the probability of an avalanche arising with sufficient charge 
to block the tubes of the quenching circuit is p then of the number Г, of ionizing parti- 
cles appearing in the counter Гр is registered immediately, and the moment at which 
the remainder L (l—p) are registered is delayed. 

After a time t = пт there still remains L(1—p) (1—p)" ~ L(1—p) e-^ unregis- 
tered discharges. If the last approximation of the equation is correct then L(1—p) — 
~ 0.65 L and р/т = A. Substituting the experimental values we have p = 0.35, 7 ~3 х 
X 107% sec, in excellent agreement with the time interval between the successive 
maxima observed on the oscilloscope (Fig. 9). Although this agreement is undou- 
btedly accidental in view of the low accuracy of the experimental evaluation of Ty this 
result points to the correctness of the estimation of the delay process for pulses in 
counters of this type. The fact that no multiple pulses as shown in Fig. 9 are observed 
to precede the pulse leading to the complete blocking of the quenching circuit valve 
can be explained by the difference between the pulses having an initial amplitude of 
several to several tens of mV and pulses blocking the circuit, the amplitude of the 
latter being of the order of 100 V in the linear region of the characteristic. Under the 
conditions in this work the simultaneous observation of the pulses of both types was 
impossible. 

From the point of view of the application of the CO, + CS, counter for geo- 
chronometric purposes, the fact that pulses occur with a delay was very undesirable 
and required additional precautions. In practice, after the counter was filled a check 
was made to see whether an increase in the length of the a. c. shielding signal above the 
standard value ry = 8.33 x 10-5 min causes a visible reduction in the number of 
pulses U. ry proved to be sufficient for all counters and filling mixtures with which 
the experiments were conducted. In fact, if the probability of an "pulse-free" dis- 
charge lasting 2 = 0.33 х 10-5 min is 0.645 SZYTA аа LOB LOA 
` then when the number of pulses in the large measuring counter oscillated by AL = 50, 
the difference in the number of background pulses resulting from the pulse delay 
effect will amount to 50 x 0.0016 — 0.08, i. e. one half of the standard error of 
a measurement lasting 50 hours at U = 100 (0(U) = 0.17). In work with small count- 
ers in Poznan the maximum difference in the counts L of a properly working appa- 
ratus was about 10. and the correction for the background increase was 0.02 as against 
a standard error of measurement of 0.09. If the measuring counter gives L cpm, then 
dL pulses arising in the, time interval between / and ¢ + dt from the moment of 
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passage of the particle will be approximately dL, = AL exp (— At) dt, and the total time 
the apparatus is blocked by these pulses will be At = tdL. Therefore the real operat- 
ing time of the apparatus is 


со 


teff = 1-- AL | «^ dt = 
0 
A: A—L : R 
and the counter efficiency 7 = 100% ° S US. With L of the order of 100 cpm (as in 
Poznan), 7 is of the order of 99.9%. 

Even for a large counter (L = 300, n = 99.5%) the inefficiency of the counter 
is of no consequence in measuring the radiation of natural carbon, but in measure- 
ments where L attains a value of several thousands, the effect can be very consider- 
able and can reduce the efficiency of the counter by several percent (Müller et al. 1950). 
In calculating the actual number of pulses from the counter, corrections for the total 
time that the apparatus was blocked by the a. c. shielding and registration signals 
was taken into account in the following way. Pulses U of a total duration С: ry, cannot 
be produced during the formation of pulses М of a total duration of Nry. It can happen, 
however, that one of the signals N occurs during the registration of an anticoincidence 
pulse in the time interval during which the apparatus is blocked. The intervals ту and 
Ty then partially or entirely coincide (Ty < Ty). The number of pulses of the а. c. 
shielding (М) falling within the total time Сту is N, = №0ту = NUty (where № 
is the number of counter pulses, N the number of signals of the arrangement). If the 
pulses appear in the time £ such that 0 <t XC Ty —Ty from the leading edge of the 
pulse U then their total time NU (ty—ty) does not participate in the blocking of the 
apparatus at all, and the rest of the pulses N, amounting to Na = N,—NU (ty—ty) = 
~NUTty block the apparatus for a mean period of ry/2. Therefore the actual counting | 
NU« м? 

2 


time for the apparatus will be 1 — Nry — Ut, + [ (ty — ty) d + 


and the real number of pulses, to a sufficient approximation, will be 
U = U[1— Uty— Ney (1 — U [ry — ty/2])] 1. 


When ry = 8.33 x 107? min and ty = 16.67 x 10-5 min (when the large counter 
was used) the mean values were М = 10? cpm, U ~ 85 cpm, so that Nty = 8.33 X 
X 107? min., Uty = 1.42 X 107? min. and the total blocking time for the apparatus 
was ~ 9.71% of the total operating time. In these conditions the effectiveness of 
counting was 90.3%. The absolute efficiency of carbon in contemporary organisms 
(17.8 F 1.2) cpm/g found in the work of Mościcki (1953) and a comparison of the total 
number of pulses of the measuring counter enclosed by material shielding filled with 
an argon — alcohol mixture (203.7 - 1.8) cpm and the number of pulses of this same 
counter filled with a mixture of CO, + CS, (209.7 F 2.2) cpm, at the same level М, 
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showed that the absolute efficiency of counters of this type in measuring the activity 
of 14C from natural carbon has a value of nearly 100%. 

c) The Influence of Contamination- by Air. 

It was tound that the contamination of the counter by air results in a lowered 
efficiency. With a leak in the counter that let in air at the rate of the order of 1 mm 


Hg/24 hrs. counts were made of L and U. (Table IV). 


Table IV 
т 

Date ae N L | U L—U, U% L, 96 
Jan. 3, 1953 592.2 72.8 12.8 59.7 100 100 
Jan. 5, 1953 598.6 63.2 TS 52.1 86.7 85.6 
Jan. 6, 1953 597.8 60.4 10.98 49.4 85.9 81.4 
Jan. 7, 1953 584.1 54.1 10.07 44.0 78.9 75.6 
Jan. 10, 1953 573.9 42.6 8.96 33.6 70.3 59.6 
Jan. 11, 1953 563.6 41.4 8.10 29.6 63.3 54.2 


Аз may be seen, the probability of the registration of particles of greater specific 
ionization and giving a constant background is larger than the probability of regis- 
tering weakly ionizing mesons. It was found that a counter showing lowered effi- 
ciency (85—90%) connected with the presence of oxygen returned to an efficiency 
of nearly 100% after a certain time and subsequently behaved like counters filled 
without any impurities. The regeneration time of counter N 37 which, as a result of 
the presence of oxygen, exhibited an efficiency of about 90% was (with the normal 
operation of the counter at about 85 cpm) 14 days, the operating time embracing 
about 150 hours. 

d. Stability of Counter Operation. 

In order to study the stability of operation of the CO, + CS, counter pulses U 
of counters NV 23 and М 37 were counted at various pressures of contemporary CO, 
from various samples and with a constant addition of about 8 mm Hg CS, and counter 
N 37 filled with carbon dioxide obtained from anthracite. The results are collected 
in Table V. It was found that the large brass counter М2 did not show any distur- 


Table V h 
Counter p/CO, mm Hg Source of Carbon ] ШЕЛ 
23 250 Contemporary 11.57-:- 0.11 10.59 = 0.34 
23 289 > 11.58 = 0.18 10.45 = 0.41 
37 278 E 11.64 = 0.11 10:03 = 0:35 
37 376.7 = | 13.03 = 0.08 11.65 = 0.51 
34 376 Anthracite 11.18 = 0.10 


*) L, was calculated by assuming 15.3 F 0.5 decays of 1*C/min. g and by taking into account 
a correction for the active volume. 


$ 
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bances in operation after counting over 5 X 107 pulses with the same filling. Counter 
М 37 gave, with different fillings, a total of over 1.5 x 107 discharges without showing 
any visible changes in its operation. No temperature effects within the range of 
0—30?C were observed. 

e. Comparison of the operating conditions of the СО, proportional counter and 
the CO, -- CS, С. M. counter working under identical conditions. 

In order to compare the operating conditions of the proportional counter filled 
with pure CO, and the С. M. counter filled with CO, + CS, a number of experi- 
ments were made in which the same counter (construction as in Fig. 2) was filled with 
pure CO, at various pressures and then with a mixture of CO, + CS, with the same 
CO, molecular pressures. A direct comparison of the operating conditions showed 
that, particularly for larger pressures (over 0.5 atm.), the proportional counters have 
a longer plateau than do the СО, + CS, С. M. counters. The thresholds of the plateau 
of the G. M. counter as a rule lie several hundred volts above the threshold of the 
plateau of the proportional counter, but below the end of the plateau of the propor- 
tional counter filled with pure CO,. Table IV gives the results obtained till now. 

According to the results obtained the possibility of avoiding the CS,, which is 
very undesirable and troublesome owing to its chemical activity, justifies the some- 
what more difficult process of purifying CO, for the proportional counter and the neces- 
sity of complicating the electrical equipment of the apparatus. 
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КРАТКОЕ СОДЕРЖАНИЕ 


b. Мостицкий, О применений счётчика Г. М. наполненного смесью СО, + CS, для 
измерения активности натурального угля. 


Часть II. Проведено ряд опытов CO счётчиками T.M., наполненными CO; + CS». 
Счётчик показывает длинное плато о малом наклоне и малой чувствительности 
на загрязнения. Изменение числа антисовпадении счётчик — экранизирующие 
счётчики в функции времени запаздывания импульсов экрана интерпретируется 
как результат запаздывания импульсов этого типа счётчиков. Согласно этой 
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интерпретации запаздывает около 65°/o импульсов CO средним временем запазды- 
вания около 0,8 m сек. Обсуждается влияние явления запаздывания на произ- 
водительность и применимость счётчика для хронометрических целей. При до- 
статочно длинном времени запаздывания импульсов экрана получается повто- 
ряющиеся результаты для предела давлений COs. Для счётчика о обьеме 2—3 л. 
получается эфективность счёта около 909/9. 
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INTENSITY RATIO IN THE *FORBIDDEN" DOUBLET 4 2958 AND 
42972 À. U. IN THE SPECTRUM OF ATOMIC OXYGEN 


Ву Г. Liszka AND H. NIEWODNICZANSKI 
Institute of Physics, Jagellonian University, Cracow 
(Received February 12, 1958) 


The aim of the reported investigations was to discover in the O I spectrum the 
appearance of the “forbidden” electric quadrupole line 4 2958 А. U. and afterwards to 
measure the ratio of the intensity of this line to that of the second line in the “forbidden” 
OI doublet, which is the previously observed magnetic dipole line 42972 А. U. 

This doublet was obtained in the radiation of a specially constructed discharge tube 
filled with a mixture of helium and oxygen. In the first part of this work optimum condi- 
tions of the discharges were found in order to obtain the much fainter component of the 
“forbidden” doublet, the line 42958 А. U. which was never observed before. In the 
second part of the work the intensity ratio of these ,,forbidden" lines in the doublet was 
experimentally found to be 1:45 -- 30%, which is in agreement in order of magnitude 
with the theoretically expected ratio 1:100. 


1. Introduction 


The fundamental electron configuration of the neutral oxygen atom is 1s? 25? 
2p*. The four equivalent 2p electrons of this configuration determine five different 
energy states, all of the same even parity, namely ЗР», *P,, ЗРо, 10, and 159, of 
which the triplet state 3P, is the ground state of the atom. According to the selec- 
tion rules no transition between these energy states is allowed for the electric 
dipole radiation. А 

Fig. 1 represents the diagram of the energy levels in О І, all belonging to the 
same ground electron configuration. Transitions which are forbidden for the electric 
dipole radiation but allowed for the magnetic dipole and the electric quadrupole radia- 
tion are shown by arrows, MD indicating the magnetic dipole and EQ — the electric 
quadrupole transitions. The figures designate the corresponding wave-lengths in 
А. U., the figures in brackets give the theoretically evaluated transition probabilities 
in sec}, The green line 4 5577 А. U. and the red triplet 4 6300, 4 6363 and 4 6292 
А. U. have been well known for many years in the spectrum of aurorae, Novae, planetary 
nebulae and night sky radiation. They were obtained much later under laboratory 
conditions. The line A en A. U. was obtained for the first time in electrical discharg- 
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es by McLennan and McLeod in 1927, the red triplet by Hopfield in 1931. The ultra- 
violet lines 2 2958 and 4 2972 A. U. in the astrophysical light sources have not been 
observed thus far, since they are situated at the border of the spectral transmittance 
of the earth's atmosphere. The line 2 2972 A. U. was obtained in laboratory condi- 
tions by Kaplan іп 1939. The line 4 2958 А. U. has been obtained for the first time 
in the present work. Hence all possible *forbidden* O I lines which correspond to 
the transition: between energy states belonging to the ground electron configura- 
tion have now been obtained in the laboratory. 


MD EQ £Q 15 


EQ EQ*MD EQ*MD 


f 2, 


2972.3|(0.090) 
2958.2 |(00008) 


(0.0078 ) 


(0.0000022 ) 
6363 


Fig. 1. Radiative transitions in the oxygen atom between energy levels belonging to the ground electron 
configuration 1522522р4, 


The experimental determination of the ratio of intensities of the components of 
the ultraviolet O I *forbidden* doublet is of great importance, since both these lines 
are due either to the pure magnetic dipole or to the pure electric quadrupole radiation 
and correspond to the transitions starting from the same initial energy level. 


In such measurements the theoretically expected ratio of the transition proba- 
bilities associated with radiations of different nature and polarity can be chequed. 


2. Apparatus and method used for obtaining the “forbidden ОТ lines 


The method applied in the present work was not a repetition of any one previously 
used for similar purposes. The apparatus was designed to enable the discharge tube 
to work as well at high current density using the electrodes, as at high frequency 
electrodeless discharges, similar to the method developed by one of us in 1933 for the 

excitation of “forbidden lines (Niewodniczański 1934). 


The discharge tube has the shape of an inverted letter 77 (see Fig. 2). The vertical 
branches of the tube are closed by glass containers enclosing cylindrical electrodes E, 
and E, made of pure aluminium foils. One side surface of each electrode has an area of 
120 cm?. In order to ensure more efficient cooling of the electrodes their outer surfaces 
are in close contact with the glass walls of containers. The horizontal part ofthe discharge 
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tube consists of a tube 45 cm long and with ап inner diameter of 0.6 cm. It has а quartz 
window Q at one end in front of which the slit of the spectrograph was placed. Owing 
to the strong heating during the discharges the whole tube was made of hard Jena 
glass and was placed in a metal vessel with flowing water. Such a shape of the discharge 
tube is convenient, since with a gas mixture of properly chosen pressure and composi- 
tion one can see in the spectrum of the discharges in the horizontal tube the spectral 
lines of atomic oxygen and helium without any stronger molecular oxygen bands. 
The dissociation of oxygen molecules takes place mainly in the vertical tubes connec- 
ting the cylindrical electrode containers with the horizontal part. Spectroscopically 


Spectrograph 


to pumps =— 


Fig. 2. General scheme of the apparatus. 


pure helium is introduced from the container Z4 in small doses through the stopcocks 
K, and K;. The part of the apparatus containing the helium is made of sodium glass 
in order to prevent the quick diffusion of helium through the walls of the container. 
Тһе discharge tube can be cut off from the pumping system by valve K,. This system 
is situated behind a liquid air trap T, and consists of a 3-stage mercury diffusion pump 
and a 2-stage rotary oil pump. The vacuum obtained was always better than 107* mm 
Hg. Oxygen was introduced from the container Z, in small doses using the stop- 
cocks K, and Kg. The container Z, is filled after opening the stopcock Кү connecting 
it with the system for oxygen production. The oxygen is obtained by heating KMnO, 
in a test tube made of resotherm glass and connected by a ground glass joint. The 
produced oxygen is transmitted through dryers 5; and 5, containing KOH and through 
the trap Т, cooled by a mixture of solid CO, and acetone. The pressure of the oxygen 
in the container was measured with a mercury manometer М. 
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At optimum conditions the discharge tube was filled with helium at a pressure 
of about 4 mm Hg and a small quantity of pure oxygen so as to maintain the ratio of 
the partial pressures of oxygen and helium equal to about 1:10. Voltage up to 6000 V 
A. C. was applied to the electrodes of the discharge tube. The power of the discharge 
current could be altered by using a variable resistor R on the primary side of the trans- 
former. The highest power applied to the discharge tube was 1.5 kW, which corres- 
ponded to a current density of 0.7 A/cm? in the horizontal tube. 


А considerable disadvantage of the apparatus was the lack of possibility of meas- 
uring the partial pressure of gases in the discharge tube. 


3. Results of the experiments 


After final adjustment of the apparatus various spectral exposures were taken. 
- The discharge tube was filled either with pure oxygen or with helium and oxygen 
mixtures of different total and partial pressures. The exposures were made using 
a high luminosity (1:5) quartz spectrograph. In spite of all precautions the spectra 
contained, besides oxygen and helium, also some Hg, Н, N and М. The source of 
these impurities was probably the KMnO, of insufficient purity that was used for 
producing the oxygen. 


The best exposures containing the sought-for line 4 2958.2 À. U. and the line 
1 2972.3 A. U. were obtained when the discharge tube was filled with helium at а 
pressure of about 3 mm Hg and oxygen at about 0.5 mm Hg and when the density 
of the current in the horizontal tube was about 0.5 A/cm?. The exposures lasted from 
1 to З hours. There was по reason to extend the time of exposures over that limit 
because after a longer discharge the “forbidden“ oxygen lines gradually disap- 
peared. This behaviour can be explained by the considerable diminishing of the 
amount of oxygen in the discharge tube owing to the oxidation of the strongly heated 
electrodes. This destroyed the optimum ratio of partial pressures of helium and oxygen. 
For this reason it would be more reasonable to apply electrodeless discharges or conti- 
nuous circulation of oxygen. 


In the close neighbourhood of the “forbidden oxygen lines which were being 
searched 6 lines were identified which were used for constructing the exact dispersion 
curve in the spectral region under consideration. These lines were: Hg — 4 2967.5 and 
4 2975.2 А. U., He — 4 2945.1 Å. U. and the second positive series of М, bands 
A 2953.2, 4 2962.0 and 4 2976.8 А. U. Since as reference lines both molecular bands 
and atomic lines were used, visual determination of their positions by means of a com- 
parator could not give satisfactory results. Therefore profiles of all reference lines as 
well as of the O I lines were measured with a microphotometer and the centres of 
gravity of all these lines were taken into consideration (see Fig. 3). The position of 
a spectral line determined in this manner contains an error of + 0.003 mm which 
gives an error of + 0.1 А for the spectral region under consideration. The positions 
of the reference lines determined in this way were used to prepare the dispersion 
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Blackening іп arbitrary units 


0.500 1.000 1.500 тт 


Fig. 3. Microphotometer records of profiles of the investigated and reference spectral lines. 


curve (Fig. 4). Hence the wave-lengths of the two “forbidden“ oxygen lines obtained 
are: 


Transition 19, — Ру: A = (2972.2 + 0.1) Å. U., 
Transition 15, — 3P,: A = (2958.1 + 0.1) А. U. 


These wave-lengths are in agreement with those evaluated from the values of terms 
given by Hopfield (1931): 4 2972.3 and 4 2958.2 А. U. 

Thus the presence in the spectrum of atomic oxygen of the *forbidden* electric 
quadrupole line 4 2958.2 А. U. has been proved. 
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Fig. 4. Dispersion curve of the spectrograph in the region of the “forbidden” ОТ lines 42972.2 and 
, 4 2958.1 A. U. 
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4. Measurements of the intensity ratio in the *forbidden* ultraviolet O I doublet 


In the experiments reported above the newly discovered electric quadrupole line 
42958 A. U was considerably fainter than the previously known magnetic dipole 
line 4 2972 А. U. Subsequent experiments were undertaken in order to measure the 
ratio of the intensities of these two lines in the *forbidden* O I doublet. 

For this purpose more than 60 spectral exposures were taken using a rhodiumised 
` seven-step filter giving in the investigated spectral region the following weakening 
ratios: 1:2.6:4.6:7.3:12.5:20.0. Since the measured lines have sufficiently close 
wave-lengths it was possible to apply the method of monochromatic photographic 
spectrophotometry. This was chequed by the parallelism of the density. curves drawn 
for both these lines. 

For the determination of the intensity ratio of these lines only 3 spectrograms 
could be used. These had an exposure time of 3 hours and a current density in the 
horizontal tube of 0.6 A/cm?, the partial pressures of helium and oxygen being 3 mm Hg 
and 0.5 mm Hg, respectively. Under these conditions the fainter line 4 2958 А. U. 
could be seen on the photographic plate in 4 steps of the filter, which made possible 
the measurement of the ratio of intensities. The mean value of the intensity ratio 
obtained in these measurements is: 


To958 : Гәдә = 1:45 


Owing to the very small intensity of the fainter component of the doublet this result 
contains a rather high statistical error of about - 30%. The order of magnitude of 
the obtained ratio is in accord with the theoretically expected value of 1:100. 
The relatively somewhat higher intensity of the electric quadrupole line may be 
caused by the specific conditions in the discharge tube. The possibility of the electric 
enhancement of some transitios in the oxygen atoms cannot be excluded. 


КРАТКОЕ СОДЕРЖАНИЕ 


Л. Лишка и Г. Неводничанский, Отношение интенсивностей в запрещенном 
дублете À 2958 u À 2972 А в спектре атомного кислорода. 


Целью работы было открытие в спектре ОГ наличия запрещенной электри- 
ческой квадрупольной линии À 2958 А и измерение отношения интенсивностей 
этой линии и другой в запрещенном дублете атомного кислорода магнитной ди- 
польной линии À 2972 А, которая наблюдалась раньше. Этот дублет был получен 
в спектре излучения специально построенной разрядной трубки, наполненной 
смесью гелия и кислорода. В первой части работы были найдены оптимальные 
условия разряда с целью получения весьма слабой составляющей запрещенного 
дублета линии ) 2958 А, которая раньше не наблюдалось. Bo второй части работы 
было экспериментально получено отношение интенсивностей этих запрещённых 


линий 1 :45+30%, которое по порядку величины сходно с теоретически предви- 
денным отношением 1 : 100. i 


"A, 
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THE APPLICATION OF 
MAGNETIC NUCLEAR RESONANCE IN FLUID HYDRODYNAMICS 


Bv A. Z. HRYNKIEWICZ 
Nuclear Physics Laboratory, Institute for Nuclear Research, Polish Academy of Sciences, Cracow 
( Received March 17, 1958) 


The recording of the nuclear resonance line by means of an mirror oscillograph 
has made it possible to employ the phenomenon of non-adiabatic transition in a stream 
of flowing fluid for studying the character of the flow. Nonadiabatic transition was used 
to mark the flowing fluid in a given section of a pipe in a time of the order of a micro- 
second. The nuclear resonance apparatus served the detector of the marked fluid after 
it has traversed a certain distance in the pipe. Іп order to illustrate the method the 
shape of pulses in the detector was studied for laminar flow and for turbulent flow in 
a smooth cylindrical pipe. 


Description of the Method 


In a previous paper (Hrynkiewicz et 41., 1957) a description was given of experi- 
ments on non-adiabatic transition in a stream of flowing fluid. The fluid, after tra- 
versing a sufficient distance in the magnetic field of an electromagnet which enabled 
it to attain a stationary value of magnetization, flowed past the poles of the electro- 
magnet into a region of a scattered magnetic field of an intensity of several tens of Oe. 
In this region a coil was wound on the pipe through which the fluid flowed (Fig. 1). A 


etectromagnet 
windings 


Fig. 1. Path of fluid flow. C — re-magnetizing coil, €. г. — resonance coil. x 
(353) 


354 l А. Z. Hrynkiewicz 

IIF o _ т ААЛА а ОЕ - -- - -- 
magnetic field directed opposite to the scattered field was set up suddenly in the coil. 
The time during which the change of direction of the magnetic field took place was 
much shorter than the period of the Larmor precession of the magnetic moments 
of the protons in this field and hence the condition was satisfied for non-adiabatic 
transition. Аз a result, the magnetization did not keep up with the change in the 
direction of the field, and in the volume of fluid subjected to sudden re-magnetization 
the magnetization became anti-parallel to the intensity of the magnetic field. The fluid 
then entered again between the poles of the electromagnet and flowed through the 
coil of the resonance apparatus. This made it possible to measure the amplitude of the 
nuclear resonance lines. The changes in the magnetic field to which the fluid was 
subjected from the moment of re-magnetization until the observation of its resonance 
lines took place slowly enough so that the condition for nonadiabatic transitions was 
satisfied. Hence the reversal of the magnetization in relation to the direction of the 
field was maintained to a considerable extent until the re-magnetized fluid entered the 
resonance coil. А drop in the amplitude was observed and at high rates of flow it was 
noted that the nuclear resonance lines reversed depending on the degree to which the 
effect was reduced by the relaxation process. The phenomenon was recorded by 
means of an mirror oscillograph. 

To the author's best knowledge, the experiment described here constitutes the 
only method of marking flowing fluid during a fraction of a microsecond in a specific 
section of the pipe without any disturbance of the flow. The resonance coil fulfils the 
role of a detector of the marked fluid and makes it possible to follow the changes in the 
shape of the marked volume of fluid during the flow. 

Ап attempt to apply the method consisted in comparison of the shape of the 
pulse recorded in the detector with the shape theoretically predicted for laminar flow 
and the illustration of the changes in the shape of the pulse of the marked fluid during 
the transition to turbulent flow after the critical value of the Reynolds number is 
exceeded. 


Theoretical Shape of Pulse for Laminar Flow 


We denote by / the distance of the front of the region in which the fluid was 
marked from the beginning of the resonance coil (Fig. 2). The volume V’ of the marked 
fluid flowing into the resonance coil is given as a function of time by 


V' = às ro ( pact 4) аг (1) 


0 


marked volume 
of liquid resonance coil 


— ДО 


44 ( 
Fig. 2. Explanation of symbols. 
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where the time 2 is counted from the moment when the resonance coil is reached by the 
front of the marked fluid flowing at a velocity of 2v, v is the mean velocity, and v the 


А 5 RAY -~ R?— r”? 
velocity of the fluid at a radius г” v = 2v m (R — radius of the pipe). After 
2 
integration and elimination of ^ by means of the relation t = E 2 we have 
29 Аа 
ERA ADU 
ЕЕ m) 


where A = xR? is the area of the pipe cross-section. Let A} be the length of the section 
of pipe in which we mark the fluid and s the length of the resonance coil. In order to 
determine the shape of the pulse in the detector (i. e. in the resonance coil) it is neces- 


=... i 


Fig. 3. Four phases of pulse recording. 


sary to find the time dependence of the volume V of the marked fluid in the coil. 
The flow of the marked fluid through the coil can be divided into four phases as illustra- 
ted in Fig. 3. 


S 
"= 
2 Av? t? 
EV = 3 
т (3) 
5 Al 
Э ске 
А 5 _ As (4% — s) 
попою» ні) обі i 
l Al 4- s 
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| Алан — Al) — (2% — $3] 
дра с (- 8. а) а 30 (5) 
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gite Li. 
PRS. к=; кау} = 2241 (6) 

The акне palet ынаны in phase Обов 
бе g- (I + AL + s? — 241] — 1) (7) 


18 
Г пах = 4 {A1 +1 5 — [0 + Al 4+ 5)? — 2015] *} 


and the maximum coefficient of filling of the degree to which the marked fluid fills 


the coil is 


V max 
"max = "ds (9) 


A plot of ава function of vt for 1 = 33 cm, ДЇ = 15 cm and s = 1.5 cm is given 
in Fig. 4. 

The derived formulae apply to the ideal case in which the pipe has the same dia- 
meter along its entire length. In the described experiments the re-magnetizing coil was 


even front case 


Pulse amplitude 


5 10 75 ст 9% 


Fig. 4. Comparison of theoretical and experimental curves. 


3 cm long but was wound onto a thickening in the pipe of an internal cross-section of 


about 5 times that of the rest of the pipe so that the volume of re-magnetized fluid. 


was equivalent to that of the above length AJ = 15 cm. In order to compare the expe- 
rimental results with the theory in this case it is better to assume that the re-magne- 
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tized fluid whose flow we are studying flows into the pipe from a reservoir. Hence the 
experimental conditions are approximated by the case in which the front of the re- 
-magnetized fluid and the front of the fluid with the normal direction of magnetization 
flowing behind it start from the same point (from the exit of the reservoir) i. e. from 
the front of the re-magnetization region at a time interval of 41/20. 

The formulae for the volume of the marked fluid for the successive phases of 
the phenomenon then take the form: 


2 Av? 12 
a) M = a i (10) 
_ As (Avt — 5) 
ы) 2-2 Q4 5) (1) 
Ср. І, Al | A[st—3 — Qx—AD 
і атавай ( 2) > | 2% + 1 саа а 


P і Al +s | Яз |4vct — s 2A1 1s — 4% 
2 чи (2 25 Е uo) 


The curve of V as a function of vt for this case and for | = 33 cm, Al = 15 cm and 
s = 1.5 cm is also represented in Fig. 4. 

The expressions given define the shape of the pulse induced in the detector by 
the marked fluid for laminar flow. Taking Re — 1160 for the critical value of the Rey- 
nolds number for a smooth cylindrical pipe, the critical velocity of the water for the 
used pipe of radius R = 0.252 cm at a temperature of 15°C is 


_ Кет 
== 


Vor = 52.5 cm/sec (14) 

‘It is interesting to compare the shape of the pulse observed experimentally for 
laminar flow with the given theoretical expressions and to show in what way the cha- 
racter of the pulse changes after transition to turbulent flow once v,, is exceeded. 


Experimental Results and Discussion 


The measurements were made for a number of mean velocities below and above 
the critical velocity. An example of recording of pulses by means of an mirror oscil- 
lograph is given in Fig. 5. For comparison of the shapes, the recorded pulses were 
normalized in such a way that their maximum amplitudes were equal. Then the curves 
were drawn as the functions of vt. For laminar flow the pulse shapes differed only 

-very slightly. In the case of turbulent flow the increase of the velocity causes the 
pulse to spread. In Fig. 4 a comparison has been made with theoretical curves of 
pulses for mean velocities of 35 cm/sec and 230 cm/sec. It is also shown the limiting 
case of the rectangular pulse that would have been recorded had the fluid flowed with 


an even front, i. e. had the flow velocity not depended on r. Experiments have shown 
ў 
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that as the velocity of turbulent flow increases the pulse spreads more and more, 
in other words, the velocity distribution in the pipe cross-section undergoes increasing 
flattening. For laminar flow, as can be seen from the figure, the shape of the pulse 
agrees quite well with that predicted theoretically. 

"The differences in the pulse rise probably come from the fact that as a result 
of the expansion and the curvature of the pipe the parabolic profile of the velocity of 
laminar flow is not formed immediately but after 
the fluid has traversed some segment of pipe. It 
would be advisable to avoid thickening and cur- 
ving of the pipe. This can be accomplished by 
employing two electromagnets. In one of them 
the fluid would be magnetized and in the other 
the resonance lines would be observed. The flow 
between the two electromagnets would be through 
straight pipe. On part of this straight linear 
section, screened to a certain degree from the 
fields of the electromagnets two coils could 
be wound, one on top of the other. One coil 
would produce a weak magnetic field parallel to 
the direction of flow and the other would serve to obtain suddenly a ma- 
gnetic field in the opposite direction. 

The drop in the experimental curve that is steeper than the theoretical curve 
is explained by the action of the relaxation process. Since the effect of the relaxation 
phenomenon on the pulse shape is difficult to grasp theoretically it would be advisable 
to employ in the experiments fluids of a longer relaxation time then water (for tap 


Fig. 5. Example of pulse recorded by 
means of an mirror oscillograph. 


water the relaxation time is T, ~ 2 sec.). 


The examples given above show that the process of non-adiabatic re-magnetiza- 
tion constitutes convenient method for marking flowing fluids, and the nuclear 
resonance apparatus can be employed for the detection of the marked fluid. From 
the shape of the recorded pulse conclusions can be drawn concerning the character 
of the flow. 

I wish to express my gratitude to Professor H. Niewodniczanki for allowing me 
to do this work in the laboratory which he heads, to Professor J. Litwiniszyn and his 
co-workers for valuable discussions, and to T. Waluga and G. Zapalski for their assis- 
tance in making the measurements. 


КРАТКОЕ СОДЕРЖАНИЕ 


А. 3. Гринкевич, Применение магнитного ядерного резонанса в динамике жидкостей. 


Регистрирование линии ядерного резонанса с помощю шлейфового осцилло- 
графа позволило применить явление неадиабатического перехода в струе теку- 
щей жидкости к исследованию характера течения. Неадиабатический переход 
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использовано для отмечения текущей жидкости B определенном отрезке трубы 
во времени порядка микросекунды. Аппаратура ядерного резонанса становит де- 
тектор меченой жидкости после пробытия нею некоторого отрезка трубы. С целью 
иллюстрации этого метода исследовано форму импульсов для ламинарного и тур- 
булентного течения в гладкой цилиндрической трубке. 
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LABORATORY EQUIPMENT AND TECHNIQUES 


A HIGH LUMINOSITY QUARTZ SPECTROGRAPH FOR THE FAR 
ULTRA-VIOLET REGION 4 2300 — 1850 À. U. 


By A. KISIEL AND H. NIEWODNICZAŃSKI 
Institute of Physics, Jagellonian University, Cracow 
( Received March 18, 1958) 


A small quartz spectrograph of great sensitivity (1:1,82) for the far ultra-violet 
22300 — A 1850 А. U. has been built at the Institute of Physics of the Jagellonian 
University. In order to avoid the influence of the absorption of atmospheric oxygen CO, 
was blown through the spectrograph. The most favourable conditions for the spectral 
fegion 42000 — 11850 А. U. are when the plateholder makes an angle of 27? with 
the optical axis of the camera lens. The mean linear dispersion at such a sitting 
is 17 А. U. per mm. 7 


Spectroscopical investigations in the far ultra-violet spectral region below 4 2300 
À. U., when using quartz spectrographs, are difficult owing to the strongly increasing 
absorption in quartz with decreasing wave-length and owing to the appearance 
of absorption bands of atmospheric oxygen (Granath 1931). The very strong absorp- 
tion of light in this spectral region in the gelatin of photographic emulsions can be 
made not effective either by the use of Schumann type photographic plates or by the 
application of the method of Ducleaux and Jeantet (1923) of sensitization of ordinary 
plates by fluorizing mineral oils. Therefore, in order to enable the use of quartz 
spectrographs in the far ultra-violet region it is necessary to fulfil the following two 
conditions. First, to make the path of light in the prism and lenses as short as possible. 
And secondly, to remove as thoroughly as possible the atmospheric oxygen from the 
entire path of light between source and photographic plate in the spectrograph. 

Rump (1925) was the first to be successful in observing the emission of the 
mercury resonance line 2 1849 AUTISM by replacing with СО, he air along the entire 
length of path of light between the resonance vessel and the photographie plate. 
According to the investigations of Schmidt (1925) carbon dioxide is transparent in 
the entire spectral region above 4 1840 À. U. A similar procedure was applied with 
success by one of us (Niewodniczanski 1928) when investigating the band absorption 
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of mercury vapour in the far ultra-violet spectral region down to 4 1849 À. U. In this 
case the stream of CO, was blown through a small quartz spectrograph, model А of C. 
Leiss, as well as through the space between its slit, the absorption vessel and the 
source of light. 

А quartz spectrograph specially designed for research work in the spectral region 
À2300— 4 1850 A. U. was built in the Institute of Physics, Jagellonian University, 
Cracow. It consists of two plano-convex quartz lenses of short focal length (соШ- 
mator lens f = 65.2 mm, camera lens f = 87.7 mm, both for A = 2310 À. U.) and 
a Cornu type double quartz prism (base length — 47 mm). 

The whole body of the spectrograph has been constructed in a form facilitating 
its filling with СО, gas, which afterwards, during the exposure, can be steadily blown 
through the spectrograph. In this spectrograph, when filled with atmospheric air the 
absorption of light at the wave length 1860 À. U. is about 76 per cent. By replacing the 
air with СО, gas the absorption is decreased to about 53 per cent, this being the ab- 
sorption in the quartz, which cannot be removed. These figures were obtained by 
using the experimental data given by Granath (1931). The introduction of CO, in 
place of air not only substantially decreases the time of exposures in the far ultra- 
-violet spectral region, but what is far more important, eliminates the necessity of 
taking into account the complex structure of the absorption bands of O, between 


4 2200 А. U. and 4 1850 А. U. 


Fig. 1. General scheme of tlie OZNA s — source, 5] — slit, D — diaphragm, Р — plateholder, 
— CO, inlet. 


The design of the spectrograph permits the angle of the plateholder to be changed 
over a vide range (20° to 45? with respect to the axis of the camera lens). The general 
scheme of the spectrograph without construction details is shown in Fig. 1. The plate- 
holder was designed for plates of the size 6.5 x 9 cm. 

The most favourable conditions for the spectrum in the region A 2000— 4 1850 
А. U. was found to be the minimum deflection angle of the prism at 4 2310 A. U. 
(incident angle of the collimator axis 53? 50") and the plateholder angle at 27°. 
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Fig. 2 represents the dispersion curve of the spectrograph corresponding to this 
adjustment. At these conditions the mean linear dispersion in the spectral region 


between 2 1860 and 21850 A. U. is 17 А. U. Йти. 
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Fig. 2. Dispersion curve of the quartz spectrograph. 
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A spectrum with well difined Al III lines 4 1854.67 and 2 1862.90 Аер 
easily obtainable with CO, -filled spectrograph on ordinary photographic plates 
sensitized with mineral oil in exposures lasting for few seconds. 


The luminosity of the spectrograph for A = 1860 A. Unis 1:182. 


КРАТКОЕ СОДЕРЖАНИЕ 


А. Кисель и Г. Неводничанский, Кварцевый спектрограф большой светосилы для 
далекой ультрафиолетовой области спектра А 2300-1850 А 


В Институте Физики Ягеллонского Университета был построен малый квар- 
цевый спектограф весьма большой светосилы (1: 1,82) для далекого ультрафио- 
лета 2300—1850 А. Для устранения поглощения света атмосферным кислородом 
спектограф продувался СО». Оптимальные условия работы в области спектра 
2000—1850 А существуют при угле наклона касеты к оптической оси объектива 
камеры равным 27°. Средняя линейная дисперсия в указанной области спектра 
при этом была равне 17 А/мм. 
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LETTERS TO THE EDITOR 


INTENSITY RATIOS OF SPECTRAL LINES IN CA I TRIPLETS 
AT DIFFERENT CONDITIONS OF EXCITATION 


By Mns. Z. LES AND H. NIEWODNICZAŃSKI 


Institute of Physics, Jagellonian University, Cracow 


( Received March 18, 1958) 


The intensity ratios of the spectral lines in the visible Cd I triplet: 


(5s5p3P9 — 5s6s3S,) 5086 AU 
(5s 5p3P9 — 55635,) 4800 AU 
(5s5p3PQ — 556535.) 4678 AU (Fig. 1) 


were investigated in a number of experimental works (Dorgelo 1924, Larche 1931, 
Ornstein 1938). Results obtained by various authors differ, however, very much from 
each other. 

In our measurements of the intensity ratios in the visible triplet in the Cd І spect- 
rum with excitation in electrodeless discharges we have found that these ratios depend 
on the conditions of excitation. In these experiments two types of sealed glass tubes 
were used: vacuum tubes with pure cadmium placed in an electric furnace for obtaining 
the required pressure of the cadmium vapour, and tubes containing cadmium in 
an atmosphere of an inert gas (argon or xenon) at a pressure of 0.1 to several 
mm Hg. In both cases for different cadmium vapour pressures different intensity 
ratios were obtained in the above-mentioned triplet, but they always deviated from the 
theoretical predictions. к 

According to the Burger-Ornstein intensity rule the intensity ratios in simple 
triplets of a sufficiently small degree of separation should be equal to the ratios of 
the statistical weights of the corresponding levels. In our case this deviated therefore be: 


ууу cou Urges LOO OO ЁЎ 
Since this triplet is considerably split the factor »* must be taken into account; hence 
we have: 

TEST ca Eo e OU 66720: 

'The measurements of the intensity ratios in Cd I visible triplet were then made by 
employing a Paschen-Schüler type hollow-cathode discharge tube as a source of light. 
In this case the measurements could be more systematic and the conditions of the 
excitation better defined/ The only variable parameters here were the current and the 
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1. Energy level diagram for cadmium atom. 


pressure of the inert gas. The cathode of the discharge tube was 40 mm in length 
and 4 mm in diameter. The emission of light took place in an atmosphere of helium 
or argon. It turned out that the intensity ratios in the given triplet change in a strictly 
defined manner with changes in both the current and the pressure of the inert gas, 
these changes being represented by a regular curve. The differences between the 
results obtained under the same conditions (the same pressure and current) are very 
slight (395). Fig. 2 shows the mean values of the intensity ratios as a function of the 
current for two substantially different pressures of helium (p, and ро). These ratios 
are so represented that the middle line 4800 À has a stable intensity of 76 and the 
intensities of the 5086 À and 4678 À lines are calculated with respect to it; according 
to the theoretical predictions they should amount to 100 and 28, respectively. As 
can be seen from the graph, the 4800 À line is always much stronger than the remaining 
two, and the intensity ratios are closer to the theoretically expected ratios for the 
smaller values of current and when the current is the same for the higher helium 
pressures. It is very difficult to measure the intensity ratios at currents below 20 mA 
because of the very small intensity of the line being studied. 

Since both states ЗРО and ЗРО on which transitions corresponding to the 5086 А 
and 4678 А lines end are metastable states, then it seems that the cause of changes 
in the intensity ratios with changes in conditions of excitation lie in the reabsorption 
phenomena connected with the changes in the population of these states. For, it should 
be noted that although both states ?P and ЗРО are metastable, they are in different 
positions with respect to} he state ?P?. Since in the conditions of our experiments 
the mean life-time of the 3P? state, as seen from simple calculations, is considerably 
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2. Intensity ratios of the spectral lines in the visible Cd I triplet plotted versus current intensity for 
two different pressures of helium in the Schüler type hollow-cathode discharge tube. 


greater than the mean time between the collisions of atoms of helium and cadmium, 

then the transition from state ЗРО to 3P? as well as from ЗРО to ЗРО could be possi- 

ble at the cost of the kinetic energy of the helium atoms. Тһе difference in energy, 

however, between the metastable state 3PQ and the state 3P? AE(3P0—3pP0) = 0.067 

eV, and the difference in energy between the metastable state ЗР? and the state ` 
зро AE(3P2—3P?) = 0.145 eV, whereas the mean thermal energy of molecules 

within the limits of temperatures 300°K — 600°К is 0.04—0.08 eV. Hence the 

exchange of energy between states ЗРО and 3P takes place continually whereas the 

probability of transition between states 3P? and 3P3 is considerably smaller. 

On the basis of these general considerations the experimental results obtained 
by the authors can be explained in the following manner: 

1. Current dependence. As the intensity of the electrical cürrent in discharges 
increases so does the number of excited Cd atoms and hence also the density of the 
absorbing layer; therefore the intensity of the spectral lines ending in metastable 
states should diminish. The intensity of the 5086 À line, as can be seen in F ig. 2, in 
fact does diminish, but as a result of the thermal exchange between energy states 
зро and 3P? an equilibrium is established between these states and the increase in 
current does not produce a change in the ratio of the populations of states 3P and 
3P0; hence the intensity ratio of lines 4678 À and 4800 À remains unchanged. 

2. Pressure dependence. Аз the pressure of the inert gas rises the exchange between 
states 3P? and ЗР becomes increasingly more probable; as a result, the number of 
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metastable states ЗРО diminishes and hence so does the reabsorption of the line 
5086 À. Therefore, with pressure increasing the intensity. of the line 5086 À with 
respect to the intensity of the line 4800 À should increase, and this does in fact happen. 
Owing to the continuous exchange of energy between the states 3PQ and ЗР? the change 
in the helium pressure does not exert any influence on the transitions between these 
states. The intensity of the 4678 À line with respect to the line 4800 À remaines cons- 
tant. 

The visible triplet Cd I was also investigated using other sources of light such as 
spark, arc and gaseous high frequency discharge tube. Different results were obtained 
for the different light sources, but they always differed from the theoretically predicted 
intensity ratios. 

All of the reported intensity measurements were made by the photographic 
method of heterochromatic spectrophotometry, employing an incandescent projection 
lamp with a known intensity distribution. 

Analogous measurements of intensity ratios were made also for the third Cd I 
triplet belonging to the same "sharp" series but lying in the ultraviolet spec- 
tral region (2868 А, 2775 А and 2734 А). The initial level in this case is 558535), 
as shown in Fig. 1. For various inert gas pressures and for different values of current 
à constant result was obtained which was in agreement, within the limits of expe- 
rimental error, with the theoretical predictions. 
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